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cn ■ Abstract 

■ Since the seminal work [30, 31] the iteration stable (STIT) tessellations have 

Q^ i attracted considerable interest in stochastic geometry as a natural and flexible yet 

I analytically tractable model for hierarchical spatial cell-splitting and crack-formation 

processes. The purpose of this paper is to describe large scale asymptotic geome- 



try of STIT tessellations in M and more generally that of non-stationary iteration 
I infinitely divisible tessellations. We study several aspects of the typical first-order 

geometry of such tessellations resorting to martingale techniques as providing a direct 
link between the typical characteristics of STIT tessellations and those of suitable 
^ \ mixtures of Poisson hyperplane tessellations. Further, we also consider second-order 

' properties of STIT and iteration infinitely divisible tessellations, such as the variance 



of the total surface area of cell boundaries inside a convex observation window. Our 
techniques, relying on martingale theory and tools from integral geometry, allow us 
to give explicit and asymptotic formulae. Based on these results, we establish a func- 
tional central limit theorem for the length/surface increment processes induced by 
STIT tessellations. We conclude a central limit theorem for total edge length/facet 
surface, with normal limit distribution in the planar case and non-normal ones in all 
higher dimensions. 
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1 Introduction 



1.1 General Introduction 

Infinite divisibility or stochastic stability of a random object under a certain operation is 
one of the most fundamental concepts in probability theory, prominent examples includ- 
ing the classical theory of infinite divisible and stable distributions with their applications 
around the central limit theorem, max-stable distributions studied in extreme value theory 
or union infinitely divisible random sets studied in the classical theory of random closed 
sets [22]. 

In the present paper we will deal with non- stationary iteration infinitely divisible random 
tessellations of the rf-dimcnsional Euclidean space and more specifically with stationary 
random tessellations that are stable under the operation of iteration - called STIT tessel- 
lations for short. In the stationary case, the purely mathematical motivation for this type 
of random tessellations goes back to R. Ambartzumian in the SOthies. The principle of 
iteration of tessellations can roughly be explained as follows: Take a random primary or 
frame tessellation and associate with each of its cells an independent copy of the tessella- 
tion itself, a component tessellation, which is also independent of the primary tessellation 
as well. Now make in each cell a local superposition of the primary tessellation and the 
associated tessellation, whereupon scale the resulting random tessellation by a factor 2 
in order to ensure that the mean surface measure of cell boundaries stays constant. The 
described operation can now be repeatedly applied and we obtain in this way a sequence 
of random tessellations. It can be shown that this sequence converges to a random limit 
tessellation and that this tessellation must be stable under iterations - a STIT tessellation 
- in the translation-invariant set-up. In the general case, the resulting tessellation is itera- 
tion infinitely divisible in any finite volume, where this property has the same relation to 
iteration stable random tessellations as infinitely divisible random variables have to stable 
ones. 

Starting with [30] , STIT tessellations and their theoretical framework were formally intro- 
duced in [31] by W. Nagel and V. Weiss. It was the same research group who discovered a 
first basic technique for studying mean values and even some distributions related to the 
geometry of STIT tessellations by writing certain balance equations based on the stochastic 
stability of the tessellation, see [32], [33] and [46] for the mean values as well as [25], [24] 
and [45] for distributional results. In [26] and [27] a new aspect was introduced into the 
theory, namely a tessellation- valued random Markov process on the positive real half-axis 
with the property that at each time the law of of the tessellation is stable under iteration. 
This process sheds light on the hierarchical and temporal structure of STIT tessellations 
and leads to a random process of cell divisions in any finite volume and moreover in the 
whole space. This point of view can be exploited to establish further results on STIT 
tessellations and in fact we will also make use of it in the present paper. In addition, 
the finite volume Markovian construction provides a link to the class of more general and 
non-stationary iteration infinitely divisible random tessellations. 

In next Subsection 1.2 we outline important facts about STIT tessellations and introduce 
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the concept of iteration infinitely divisible random tessellations. Afterwards, in Subsection 
1.3 the plan of the paper as well as a survey of our main results are presented. 

1.2 Iteration Infinitely Divisible and STIT Tessellations 

A tessellation of M*^ is a locally finite partition of the space into compact convex polytopes, 
the cells of the tessellation. One can regard a tessellation either as a collection of its cells or 
as the closed set of their boundaries. We will mostly follow the second mentioned path and 
denote by Cells (F) the set of cells of the tessellation Y (by the Jordan-Schonflies theorem 
the correspondence between Y and Cells(F) is one-to-one). Thus, a random tessellation 
can be regarded as a special random closed set in the sense of [40]. In particular, this 
imposes the usual Fell topology and the corresponding Borel measurable structure on the 
family of tessellations, see ibidem. A random tessellation Y (regarded as a random closed 
set in Mf^) is called stationary if its distribution does not change upon actions of transla- 
tions. Analogously a random tessellation is called isotropic if its distribution is invariant 
under the action oi SO{d). 

Whenever wc have two random tessellations Yi and Y2 of we can define their itera- 
tion/nesting. To do so, we associate to each cell c G Cells(Yi) an independent version ^2(0) 
of Y2 and we assume furthermore the family {Y2{c) : c e Cells(yi)} to be independent of 
Yi. Then we can define the iteration of Yi with Y2 by 

yiffly2:=l'iffl{i^2(c) :ceCells(yi)}:=yiU U (Y2{c)r\c), 

ceCeiis(yi) 

i.e. we take the local superposition of Y2 and the family {12(c) : c e Cells(yi)} inside the 
cells of Yi. It was shown in [27] that Yi SY2 is a stationary random tessellation as soon 
Yi and Y2 are. A stationary random tessellation Y is called stable under iteration or STIT 
for short iff 

m (;y ffl. . .ffly)^ ^y, m = 2,3,... , (i) 

m 

where = stands for equality in distribution, i.e. if its distribution does not change under 
rescaled iteration. In fact, using the uniqueness results, see Theorem 3 and Corollary 2 in 
[31], it is easy to see that it is enough to take one fixed m > 1 in (1). 

To proceed, let us be given a constant < t < 00 and a probability measure TZ on the unit 
sphere Sd-i usually identified with the induced distribution of orthogonal hyperplanes on 
the space "Ho of (d — l)-dimensional linear hyperplanes in M.'^, also denoted by TZ in the 
sequel for notational simplicity. Define the measure A on the space T-L of affine hyperplanes 
in as the product measure 

A:=i+®n (2) 
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Figure 1: Realizations of a planar and a spatial stationary and isotropic STIT tessellation 
(kindly provided by Joachim Ohser and Claudia Redenbach) 



of standing for the Lebesgue measure on the positive real half-axis (0, oo), and of 7^, 
where a pair (r, m) G (0, oo) x Sd-i is identified with the hyperplane {x G M*^, {x,u) = r}. 
Throughout this paper we always require that the support of 71 spans the whole space, i.e. 

span(supp(7^)) = M"^. 

Assume now that we are given a stationary random tessellation Y with surface intensity 
t (i.e. the mean surface area of cell boundaries per unit volume equals t) and directional 
distribution TZ (i.e. the distribution of the normal direction of the face containing the 
typical point is given by 71) and define the sequence (X„(F)) by 

Ii(Y) := 2{Y ffl Y), ln{Y) := X„_i(F) ffl nF = n (F ffl . . . ffl F), n > 2. 

n — 1 ^ ' 

n 

It was shown in [31, Thm 3] that X„(y) converges in law, as n — )■ oo, to a stationary 
random limit tessellation Y{tA) uniquely determined only by tA. This tessellation is easily 
shown to be stable under iterations and is called the STIT tessellation with parameters t 
and A. Without confusion we will write Y{t) instead of Y{th), whenever the measure A or 
7Z is fixed. It is a crucial feature of Y{t) that it admits a very natural and intuitive explicit 
construction. For a restriction Y{t,W) of Y[t) to a compact convex window W G 
this construction can be informally be described as follows (the reader is referred to [31] 
for full details). Assign to the window W an exponentially distributed random lifetime 
with parameter A([iy]) where \W] := {H e Ti, H H W 0} stands for the family of 
all hyperplanes hitting W. Upon expiry of its lifetime, the cell W dies and splits into 
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two sub-cells W'^ and W~ separated by a hyperplane in [1^] chosen according to the 
law A(-)/A([W^]). The resulting new cells W'^ and W~ are again assigned independent 
exponential lifetimes with respective parameters A([W^"'"]) and A([W^^]) (whence smaller 
cells live stochastically longer) and the entire construction continues recursively, until the 
deterministic time threshold t is reached. The cell-separating (d — 1) -dimensional facets 
arising in subsequent splits are usually referred to as I-faces (or I-segments for d = 2 as 
assuming shapes similar to the letter I). The described process of recursive cell divisions 
is called the Mecke-Nagel- Weiss- or (MNW)-construction in the sequel and the resulting 
random tessellation created inside W is denoted by Y{t, W) as mentioned above, whereas 
the collection of all I-facets or I-segments is denoted by MaxFacets(F(t, 1^)). Moreover, 
we write MaxFaceSfc(y(i, 1^)) for the collection of /c-dimensional I-faces of Y{t, W), where 
by a /c-dimensional 1-face we mean the maximal union of connected and fc-coplanar k- 
dimensional faces. In fact, A:- dimensional I-faces ofY(t, W) can also be alternatively defined 
as the A;-faces of I-facets. It was shown in [31] that the law of Y(t, W) is consistent in W, 

i.e. Y{t, W)r\V = Y{t, V) for convex V <Z W and thus Y{t, W) can be extended to random 
tessellation Y{t) — Y{tA) on the whole space, which is then proved (cf. [31]) to coincide 
with the hmit tessellation Y(t) considered above, as the notation suggests. Again, the sets 
of all I-faccts and I-faccs of Y{t) arc denoted by MaxFacets(y(t)) and MaxFacesfc(y(t)) 
(0 < k < d — 2), respectively. The stationary random tessellation Y{t) is additionally 
isotropic if and only if TZ is the uniform distribution on Sd-i- In this case we will write 
Ajso or A^Q for the invariant measure on the space of hyperplanes H. 
A simple yet crucial observation is that even though only translation-invariant measures A 
of the form (2) show up in the limiting STIT tessellations, the MNW-construction can be 
carried out with arbitrary non-atomic and locally finite driving measure A (i.e. A([iy]) < oo 
for W bounded) on "H also leading to a consistent family Y{t, W) and eventually, by 
extension, yielding Y{t). Many of our theorems below will be stated in this general context. 
It should be emphasised though that such tessellations are no more iteration stable (STIT). 
However, they have the general property of being iteration infinitely divisible, as they can 
be readily checked to arise as m-fold iterations of Y{t/m) for each m > 2 in all finite 
volumes. Formally, this means that 

Y{t/m,W)^'^ ^Y{t,W) 

for all compact convex sets 1^ C R*^, which follows directly by the MNW-construction as 
yielding 

y(s, W) ffl Y{u, W) = Y{s + u,W). 

It is more than natural to expect that also Y{t) = Y{t/m)^"^ for all m in the whole 
which should be easily provable by adopting the theory developed in [27], thus even better 
justifying the term iteration infinitely divisible in our context, yet this falls beyond the 
scope of the present work. 

The STIT tessellations enjoy a number of remarkable properties. Below, we recall some of 
them which will be of importance for our further argumentation. All proofs can be found 
in [30] and [31]. 
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(a) STIT tessellations have Poisson typical cells, i.e. under stationary A the interior of 
the typical cell of Y{t) - denoted by TypicalCell(y(t)) - has the same distribution as 
the interior of the typical cell of a stationary Poisson hyperplane (or line or plane) 
tessellation with the same surface intensity t > and the same directional distribution 
TZ, see Lemma 3 in [31]. The difference between these two types of tessellations arises 
from the mutual arrangement of the cells, see Fig. 1. 

(b) Sectional STIT tessellations are STIT tessellations as well. In particular, let 

be a /c-dimensional plane in R'', 1 < k < d — 1, and let Y{t) be a stationary and 
isotropic random STIT tessellation with parameter t > 0. Then Y{t) fl has the 
same distribution as the stationary and isotropic STIT tessellation with parameter 
Xkt constructed inside the plane E^, i.e. 



Y{tA^,,R'^)nE,^Y(XktAil,E,) 

with, cf. [28, (3.29T)] 

by abusing the above introduced notation, where we have assumed the window to 
be compact. Moreover, if A; = 1 the intersection of Ei with Y{t) induces a motion- 
invariant Poisson process with intensity 

r (^) 



V^T (^) 

on the line Ei. In view of the above considerations we see that the motion- invariant 
STIT tessellations in dimension 1 are the motion-invariant Poisson processes (note 
that in dimension 1 motion-invariance is the same as translation-invariance). Thus, 
STIT tessellations can be seen as another generahzation of the 1-dimensional Poisson 
process to higher dimensions beside other models, see [29]. 

(c) STIT tessellations have the following scaling property: 

tY{t)^Y{l), 

i.e. the tessellation Y{t) of surface intensity t upon rescaling by factor t has the same 
distribution as Y{1), the STIT tessellation with surface intensity 1. 

We will use in the paper the same notation for iteration infinitely divisible random tes- 
sellation and STIT tessellations hopefully without confusion, mentioning in any case if we 
are working with a general locally finite non-atomic measure A, a translation-invariant one 
or even with Aj^o- 

We will also make use of the following notation in the paper: 
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• Br — B'^{o) is the d-dimensional ball around the origin with radius > 0. 

• The /c-dimcnsional volume measure will be denoted by Volfe. 

• Kj :— Wo\j{B{) is the volume of the j-dimensional unit ball, juj its surface area. 

• The uniform probability measure on the unit sphere S^-i in (normalized spherical 
surface measure) is denoted by 

1.3 Short Overview of the Paper and the Results 

As already remarked above, a STIT random tessellation induces a continuous time tessel- 
lation-valued Markov process on the positive real half-axis. For such processes, the notion 
of a generator is available. Using standard theory of Markov processes, several martingales 
associated with the tessellations under consideration can be constructed. This fundamen- 
tal observation will be exploited in Section 2 and most of our results will be based on it. 
In Section 3 we mainly deal with first-order properties and the typical geometry of iteration 
infinitely divisible or more specifically stationary and stationary and isotropic STIT tessel- 
lations. We generalize some properties known for STIT tessellations to the non- stationary 
case and derive explicit formulas for mean values and higher moments of /c-dimensional 
I- faces in the isotropic case (Sections 3.2 and 3.3). In particular we show that the distri- 
bution of the (time-marked) typical /c-dimensional I-face is a mixture of suitably (marked 
and) rcscaled /c-dimensional Poisson cells. As a main tool we use the martingale techniques 
introduced before and we compare the tessellations with certain mixtures of Poisson hyper- 
plane tessellations (Section 3.1). Typical geometry beyond the isotropic regime is explored 
in Section 3.4, where for example the conditional distribution of the typical A;-dimensional 
I-segment given its birth time for all space dimensions is determined and some global and 
local mean value formulas for the stationary as well as for the non-stationary case are 
derived. 

Second-order parameters of iteration infinitely divisible and STIT tessellations are the con- 
tents of Section 4. Based on a speciahzation of our martingale technique (Section 4.1), we 
calculate the variance of a general face-functional and as a special case we find the variance 
of the total surface area of cell boundaries in a bounded convex window. The resulting 
integral expression can be explicitly evaluated in the stationary and isotropic case by ap- 
plying the affine Blaschke-Petkantschin formula (Section 4.2). For the particular case of 
space dimension 3, an exact formula without further integrals is presented as well as general 
variance asymptotics relying on techniques from integral geometry. Certain chord-power 
integrals will reflect the influence of the geometry of the observation window to this vari- 
ance asymptotics. The purpose of Section 4.3 is to compare our martingale approach with 
a recent elegant second-order theory for stationary and isotropic planar STIT tessellations 
independently developed by W. Weiss, J. Ohser and W. Nagel. As a by-product, we derive 
there an explicit expression for the pair-correlation function of the random surface area 
measure in general space dimensions. For space dimensions 2 and 3 we also calculate the 
so-called iC-function. 
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Following the afore-mentioned first- and second-order theory we turn in Section 5 to the 
central limit problem for STIT tessellations. In Section 5.1 we show the convergence of 
the surface/length increment process wrt. to some initial time instant associated with 
STIT tessellation to a time changed Wiener process, the proof of which is based on the 
ergodicity property of STIT tessellations and the functional central limit theorem for 
square-integrable martingales. It is interesting to see that in terms of the time incre- 
mental MNW-construction of STIT tessellations our central limit theorem is more reliant 
on independencies arising in construction time than those of spatial nature. Even more 
interestingly, the 'big-bang' phase near time zero turns out to play a crucial role. In Sec- 
tion 5.2 we estabhsh normal convergence of the total edge length in the planar case and 
in Section 5.3 it is shown that non-normal limit distributions appear in the limit for space 
dimensions greater or equal 3. The argument goes by using our variance calculations from 
Section 4 and the observation that for space dimensions > 3 the initial 'big-bang' phase of 
the MNW-construction brings a non-negligible contribution to the variance of the surface 
increment process, whereas in dimension 2 this contribution is negligible but already with 
an extremely slow rate of decay. 

The aim of the last Section (Section 6) is to put our results in a more general context 
by comparing them with the results for Poisson hyperplane tessellations and for Poisson- 
Voronoi tessellations with the same surface intensity. 

In all cases, our general statements are illustrated by examples in space dimensions d — 2 
and d = 3, sometimes leading to already known formulas but often extending them at least 
in some aspects. 

Our results could - beside their intrinsic mathematical motivation - be of potential interest 
for applications for example in context of questions concerning statistical model fitting of 
random tessellations to real data. Central hmit theorems could be a basis for statistical 
inference of tessellation models and related functionals. Asymptotic confidence intervals 
and tests with respect to mean values can be derived from them, since we can make the 
first- and second-order moments explicitly available. 

2 Martingales Associated with Iteration Infinitely Di- 
visible or STIT Tessellations 

The finite volume continuous time incremental MNW-construction of iteration infinitely 
divisible random tessellations or more specially stationary STIT tessellations, as discussed 
in Section 1.2 above, clearly enjoys the Markov property in time parameter, whence natural 
martingales arise which will be of crucial importance for our further considerations. To 
discuss these processes we need some additional terminology and notational conventions. 
First, we fix a compact convex window W and a general diffuse (non- atomic) and locally 
finite measure A on %. Next, whenever a new facet / lying on a hyperplane H is born in 
the course of the Mecke-Nagel- Weiss- (MNW)-construction, that is to say a cell splits, we 
declare that its initial/birth point is the point (vertex) with lowermost first coordinate, to 
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be denoted in the sequel (in fact any other deterministic and measurable algorithm 
of distinguishing one point of / would do as well in this context). Moreover, for a given 
tessellation Y (usually taken to be Y{t, W) for some t > 0) we also introduce the notation 
i{H; Y) for H - a hyperplane in [W], with the following meaning 

• The hyperplane H is tessellated by the intersection with Y, and thus split into a 
number of {d — l)-dimensional polyhedra /i, /2, • • • which could potentially become 
new facets for the tessellation, were its construction continued. 



Whenever x G l{H;Y), we define Facet{x, H\Y) to be the new facet to be added to Y 
should a facet birth (cell split) occur at a; on if in Y. Moreover, we denote by Cell(a;, H\Y) 
the cell split by the birth of Facet(x, H\Y) and we write Cell+(a;, H\Y) and Ceir(x, H\Y) 
for the two sub-cells into which Cell(x, H\Y) divides, lying on the positive and negative 
sides of H, respectively. 

With this notation, it is easily seen that for a fixed measure A, {Y{t, W))t>o is a pure jump 
Markov process with values in the space of tessellations of W and with the property that 



dF(Y{t + dt) = YUF&cet{x;H\Y)\Y{t) = Y) = l[x G L{H;Y)]A{dH)dt, (4) 
with H e[W], and 



with \i{H; Y) \ standing for the cardinality of l{H; Y). Indeed, this is because, conditionally 
on Y{t) — y, during the period (t, t + dt] of the MNW-construction we have, as discussed 
in Subsection 1.2, 

• For each cell c G Cells(y) the probabihty that it undergoes a split is A([c])di, 
moreover the probability that two or more cells split is o{dt), whence the proba- 
bility that no split occurs is 1 — ^X]ceCeiis(y) -'^(H)) ^ (^)- Observing that 
\l.{H; Y)\ coincides with the number of cells in Y that H intersects, we get in addi- 
tion X^ceCeiis(y) -^([^D ~ I\w] which yields the remaining equality in 



• Should a cell c G Cells(y) split, the splitting hyperplane H is chosen according to the 
law 1[H G [c]](-)/A([c]), whence the probability of observing a split of c induced by 
H during the time period (t,t + dt] is just 1[H & [c\]{dH)dt. Now, having x G l.{H; Y) 
is by definition equivalent to there being a cell c = c{x, H) G Cells (F) with H G [c] 
and x — i{c f] H). For so chosen c we see that on the event Y{t) — Y, having 
Y{t + dt) = Y U Facet(a;; H\Y) is equivalent to having c split by H during (t, t + dt\. 
As noted above, the latter happens with probabihty A{dH)dt whence (4) follows. 



. We put L{H-Y):={t{f,),t{f2),...}. 





(5). 
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It should be emphasised at this point that the differential notation in (4) and (5) employing 
the symbols d¥ and dt is widely accepted in the theory of pure jump continuous time 
Markov processes and makes perfect formal sense as a commonly recognised abbreviation 
for the usual description in terms of waiting times, with dF{Y{t + dt) = Y'\Y{t) = Y) = 
a{Y, Y')dt understood as ^ while in state Y, wait an exponential time with parameter a{Y, Y') 
and then jump to Y' unless some other jump has occurred prior to that . We refer the 
reader to Chapter 15 in [5] and especially to Section 15.6 there, where this construction 
is formalised with full mathematical rigour in terms of waiting times as noted above. As 
readily verified, specialising the generic waiting-time construction to the case of (4)) and 
(5) yields precisely the standard MNW construction. 

Using (4) and (5) we conclude by general theory of Markov processes and their infinitesimal 
generators, see Chapter 1 in [20] or Chapter 15 and especially Sections 15.4 (Def. 15.21) and 
15.6 in [5] specialised for the pure jump case, the generator for {Y{t, M^))t>o is L '■— "^K-yv 

hF{Y)= [ V [F{YUFeicet{x;H\Y))-F{Y)]A{dH) (6) 

for all F bounded and measurable on space of tessellations of W. Consequently, again by 
standard theory as given in Lemma 5.1 Appendix 1 Sec. 5 in [18], see also Section 1.5 in 
[20], or alternatively by a direct check straightforward in the present set-up, we readily see 
that for F bounded and measurable the stochastic process 

F{Y{t,W))- [ hF{Y{s,W))ds (7) 
Jo 

is a martingale with respect to the filtration generated by {Y{s, Vr))o<s<t- More gener- 
ally, for bounded measurable G = G{Y, t), considering the time- augmented Markov process 
{Y{t,W),t)t>o and applying usual theory, see again Lemma 5.1 in Appendix 1 Sec. 5 in 
[18], or simply by performing a direct check, we see that 



G{Y{t,W),t) - £ (^[hG{;s)]{Y{s,W)) + ^^GiY(s,W),s)^ ds 



(8) 



is also a martingale with respect to as soon as G{Y, t) is twice continuously differentiable 
in t and supy^ \^G{Y,t)\ + \^G{Y,t)\ < +00, which is condition (5.1) in [18, App. 1 Sec. 
5]. 

To proceed, consider F of the form 

E,(y):= Yl '^(/) (9) 

/eMaxFacets(y) 

where, recall, MaxFacets(y) are the maximal facets of Y (the I-segments in the two- 
dimensional case) whereas </){■) is a generic bounded and measurable functional on [d— 1)- 
dimensional facets in W, that is to say a bounded and measurable function on the space 
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of closed {d — l)-dimensional polj^opes in W, possibly chopped off by the boundary of W, 
with the standard measurable structure inherited from space of closed sets in W. Whereas 
the so-defined F is not bounded and thus (7) cannot be applied directly, we can apply it 
for Fm := (F A iV) V -N, N e N which is bounded and let iV ^ oo to conclude that 
F{Y{t,W)) - J^hF{Y{s,W))ds is a local $>t-niartingale, see Definition 5.15 in [19] and 
take Tn — inft>o \F(Y(t, W))\ > N there. Now, apply the proof of Lemma 1 in [31] where 
the number of cells in Y{t, W), and hence for all Y{s, W), s < t, is bounded by a Purry- 
Yule-type linear birth process whose cardinality at any given finite time admits moments 
of aU orders, to conclude that {F{Y{t, W)) - LF(F(s, W))ds)t<a is of class DL for all 
a > in the sense of Definition 4.8 in [19]. Using now Problem 5.19 (i) in [19] we finally 
conclude that F{Y{t, W)) — L,F{Y{s, W))ds is a martingale. Thus, applying (6) and (7) 
for F = we see that 

E^{Y{t,W))- [ [ J2 <P{Facet{x, H\Y{s,W)))A{dH)ds (10) 
is a martingale with respect to ^jj. 

3 First-Order Properties and Typical I-Faces 

In this section we establish a number of first-order properties of Y(t,W) — Y{t,A,W) 
for general locally finite non-atomic measure A, essentially obtained by comparison with 
suitable mixtures of Poisson hyperplane tessellations. Many of these properties are known 
in the translation-invariant set-up. Afterwards, these results are used to calculate several 
mean values for stationary and stationary and isotropic STIT tessellations Y{t). 

3.1 Distribution of Cells and Lower Dimensional Faces 

The key to our results is formula (10) from Section 2. To exploit it, consider the random 
measures 

^/V^nW:^ J2 4, W := EA4^(*'^) (11) 
ceCeiis(y(t,w)) 

with Sc standing for the unit Dirac mass at c. In full analogy, define M.^^'^^*^'^^ and 
iyjPHT(tA,w) ^jigpg PHT(iA, W) is the Poisson hyperplane tessellation with intensity mea- 
sure tA, restricted to W. Further, put 

/eMaxFacesfc(y(t,W)) 

where, recall, MaxFaceSfc(y) is the collection of /c-dimensional I-faces of Y and, hkewise, 
define 

j:PHTitA,w) J2 5f, Ff := EJ-f'«^(*^'^), k^l,...,d-l. 

/eFacesfc(PHT(tA,W)) 
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Our first claim is 
Theorem 1 We have 

M^(t,W) ^ ]^PHT(tA,W)_ 

Note, that in the particular case of A being translation-invariant this reduces to the known 
fact that the typical cell of Y{t) = Y{tA) coincides with the typical cell of PHT(tA). 

Proof of Theorem 1 Using (6) and (7) with 

F{Y):= Yl 

ceCeiis(y) 

for general bounded measurable cell functional (f), with localization argument as the one 
preceding (10) we conclude that 

/ - r / Yl [<P{Cell^{x, H\Y{t,W))) + 

J Jo J[W]^^^^jj;Y) 



(j){Cen.-{x, H\Y{t, W))) - 0(Cell(a;, H\Y{t, W)))]A{dH)ds (13) 

is a Q't-martingale. For a polyhedral cell c C W. possibly chopped off by the boundary of 
W, and for H e [c] we write c^{H) and c~{H) to denote the cells into which c gets divided 
by H, lying respectively on the positive and negative side of H. With this notation, (13) 
says that 

/ f [ [ [(j){c+{H)) + (l){c-{H)) - (l){c)]A{dH)M''^''^\dc)ds 

J Jo J J[c] 

is a ^Jj-martingale. Taking expectations leads to 

/" 0dM^(*'^) = f [ [ [(j){c+{H)) + (j){c-{H))-(j){c)]A{dH)M''^''^\dc)ds 

J Jo J J[c] (•3^4-) 

for all bounded measurable 0. To proceed, we regard M^*^*'*^) as an element of the space of 
bounded variation Borel measures on the family of polyhedral sub-cells of W endowed with 
the standard measurable structure inherited from the space of closed sets in W. Consider 
the linear operator T\ on this measure space, given by 

7a(/x)= / / [5,+^H) + Sc-iH)-Sc]A{dH)pi{dc). (15) 

J J[c] 

By the definition (15), \\TM\\tv < (/[w] II Hiw]) II 

^IItv where ||-||rpY is the 

standard total variation norm of a measure, see [3, Def. 3.1.4]. This inequality turns 
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into equality when /i — 5w- Consequently, T\ is a bounded operator of operator norm 
A([l^]) < +00. The relation (14) can be rewritten in differential form 

^M^it,w) ^ TaM^(*'^'), M^(0'^) = Sw (16) 

CLv 

which, in view of the above properties of T\, admits by standard theory (cf. [17, IX. §2, 
Sec. 2]) the unique solution 

_ exp{tT,^)dw, t > 0. (17) 

It is easily seen that exactly the same equations (14), (16) and thus also (17) hold for 
I^PHT(tA,w)_ particular, M^(*'^) = M^^^^*^-^) as required. This completes the proof. □ 

Having characterized M^*^*'^) we now turn to F^*^*'^^. 
Theorem 2 For all k — 1, . . . ,d — 1 we have 

t 



Jo s 



Proof of Theorem 2 Fix /ce{l,...,d — 1}. Let -0 be a general bounded measurable 
function of a /c-dimensional I-face, as usual regarded as a closed subset of W, and for a 
{d — l)-dimensional I-facet h put 

m--- E ^(/) (18) 

/eFacesfc(/i) 

noting that the /c-dimensional I-faces of the tessellation Y{t, W) are precisely the /c-faces 
of its I-facets. Using (10), taking expectations and recalling (11) we see that 

EE^(y(t,iy)) = / (Pd¥^[\^^ = f ! ! <P{cn H)A{dH)M''^''^\dc)ds. 

J Jo J J[c] 

Applying Theorem 1 we get 

f(pd¥^lT^= f [ [ (j){cnH)A{dH)Wf''^^'^''^\dc)ds. 

J Jo J J[c] 

However, applying Slivnyak's theory, see e.g. [41, Thm 1.15], we obtain 

(/.(cn//)A(di/)MP«T(^'^)(dc)= / 0(iFf;'(^'^) 



and thus, upon taking sA in place of A, more generally, 

J ^ 0(cnif)A(dif)MP"T(^^'^)(rfc) = -^J 0rfF™^(^^'^) 

13 



whence 

follows. Now note that, by (18), 

because each k-dimensional I-face is a k-face of precisely one I-facet in Y{t, W). Moreover, 

because each k-facc of PHT(sA, W) is a k-face of {d - k)2'^-^-^ facets of PHT(sA, W), see 
Theorems 10.1.2 and 10.3.1 in [40]. Hence, we conclude that 

j = {d- k)2'^-'^-' j'^-J V^dFf^^^^'^^d. 

for all bounded and measurable, which completes the proof of the Theorem. □ 

Some of our argument in the sequel will require a straightforward formal extension of 
Theorem 2. Namely, we formally mark all 1-faccts of the tessellation F (t, W) by their hirth 
times. This gives rise to the birth-time augmented tessellation Yit^^W) with birth-time- 
marked I-facets and makes the MNW-construction of y(t, W) into a Markov process whose 
generator L is a clear modification of L as given in (6): 

LF(y) = / V [F{Y U [Facet(x; H\Y), s]) - F{Y)]A{dH) 

J\w] Jtt!^.^ 



xet{H;Y) 



for F bounded measurable on the space of birth time-marked tessellations of W. Conse- 
quently, writing F^^*'^'*, A; = l,...,d— 1, for the birth-time-marked version of F^*^*'^^ 
where each k-dimensional 1-face is marked with its birth time, by a straightforward modi- 
fication of the proof of Theorem 2 we are led to 



Corolleiry 1 For all k — 1, . . . ,d — 1 we have 



Fl^'^"^^ ^ {d - k)2'-'-' 



ds. 



3.2 Distribution of typical I-Faces and Related Mean Values 

We are now going to apply the results obtained in the last section to the stationary and 
isotropic set-up, i.e. with A = Aj^o and to calculate the distribution and the mean /c- volume 
as well as other related mean values of the typical /c-dimensional I-face ijf^ , k — 0, . . . , d—1, 
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Figure 2: Vertices (left hand side) and crossing nodes (right hand side) that appear in 
3-dimensional STIT tessellations 



of the STIT tessellation Y{t) in W''. In this paper we will restrict our attention only to 
proper A;-dimensional I-faces, which we call k-dimensional I-faces for short, see Section 
1.2 above for definitions. However, it should be noted that starting with dimension 3, 
there appear also the so-called crossing k-faces arising as intersections of traces of higher- 
dimensional I-faces inside another face. For d = 3 these crossing k-faces appear only in 
the form of nodes. The 0-dimensional I-faces are simply the vertices and the other type of 

0- faces (nodes) are called crossing nodes. Figure 2 should clarify the situation. In higher 
dimensions the structure of crossing k-faces is even more complicated, with the notable 
exception of the case k = d — 1 where there only exists one type of I-faces, namely the 

1- facets. As stated above, throughout this paper whenever talking about /c-dimensional 
I-faces we always mean the proper ones. To proceed, let (pk '■ MaxFaces^ — )■ M be a 
translation- invariant, non- negative measurable function ior 1 < k < d — 1 and denote by 
'(p^{Y{t)) the (possibly infinite) (/7fe-density of Y{t) in the sense of [40, Chap. 4.1], i.e. 

^ /eMaxFaceSfe(y(t,rVF)) 

with C some bounded convex set with positive and finite volume. The existence 
of this limit is guaranteed by Thm 4.1.3 ibidem. Using now Theorem 2 from above and 
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Campbell's formula, we obtain, possibly with both sides infinite, 

^ ^ /eMaxFacesfe(F(t,rW)) 



r-).oo r"'Vold(VV) Jo U 



lim — — -E V u>k(f) 

^ ' /eFacesfc(PHT{s,rW)) 



ds 



= id- fc)2'^-'=-^ / - 
Jo s 

= {d-k)2''-'-' f-^^{VWI{s))ds, (19) 
Jo * 

where PHT(s) (PHT(s, rW)) is the stationary and isotropic Poisson hyperplane tessellation 
with surface intensity s (restricted to the window rW). We use the formula (19) with (/jj. = 
1 and denote in this case Tp^{Y{t)) =: Nfl, referred to as the intensity of k- dimensional 
I-faces. Then we have upon applying [40, Thm. 10.3.3] together with (19) 

<] = (-^)2-/'i(:)-(t^)^^<'-«'-'=)^"-'?0(9^)^;„, 

For a /c-dimcnsional polytope /, let c(/) be some associated center function with the 
property that c(/ + x) = c{f) + x for any x & (take for example the Steiner point or 
the center of the minimal circumscribed ball of /). Now, we use again (19) with (pk{f) '■— 
![•](/ - c(/)) together with [40, Eq. (4.8,4.9)] and obtain in this case 

<]Qr<'' = (d - ^2^-'-' /' i7ST,.,«r"<"'*'-. (21) 

JO * 

where 7pht(s) mean number of A;- faces of PHT(s) per unit volume, Q|[^"^^*^ is the 

distribution of the typical /c-face of PHT(s) and Q^^*^ is the distribution of the typical 
A;-dimensional I-face of the STIT tessellation Y{t), see [40, Thm. 4.1.1] and the discussion 
ibidem for definition of typical faces and typical grains of particle processes in general. 
Inserting now the value for Njf'j and the one for 7pht(s) from [40, Thm. 10.3.3] and noting 

that Nj^'^j — ^2'^~*^~^7^^j^rp(j) by comparing the total masses in both sides of (21), we arrive 

at 

Theorem 3 The distribution Q^^*'' of the typical k- dimensional I-face ofY{t) is given by 



Y{t) /■* ds'^ ^^PHT(s) , 
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where Q^. is the distribution of the typical k-face of the stationary and isotropic Poisson 
hyperplane tessellation PHT(s) with surface intensity s. 

The preceding theorem can be rephrased by saying that the distribution of the typical 
k- dimensional I-face of a STIT tessellation is a mixture of suitable rescalings of the distri- 
bution of the typical fc-dimensional face of Poisson hyperplane tessellations. 
In what follows, ijf^ stands for the fc-polytopc with distribution Q^''*^. For the typical 
I-segment of the stationary and isotropic STIT tessellation Y(t) in M*^ we can obtain 
even more than the statement of Theorem 3, namely the fact that the length distribution 
of /j'^'' is a mixture of exponential distributions, answering thereby a question formulated in 
the Outlook of [24]. It is also interesting to note that formally marking the A;-dimensional 
I-faces with their birth-times and repeating the argument leading to Theorem 3 with Theo- 
rem 2 replaced by its time-marked extension in Corollary 1 we obtain the birth time-marked 
extension of Theorem 3 

^ Y(t) 

Corollary 2 The distribution Qj. ^/ typical birth-time-marked k-dimensional I-face 
ofY{t) is given by 



Ant) l-ds'-'^ 



^0 



^PHT(.) 



ds. 



- ' t^ 

To proceed, using Theorem 3, for the length density pf^{x) of the typical I-segment of a 
stationary and isotropic STIT tessellation Y{t) in we obtain because of key property 
(b) from Section 1.2 

pf\x) = f X.se^'^^^^ds = jT-^^l{d+l,Xitx), (22) 







td (Xityxd+^ 



where 7(-, •) is the lower incomplete Gamma-function and Ai is given by (3). Prom this it 

is easily seen that for the length of the typical I-segment only the moments of order 1 to 
d — 1 exist, see Section 3.3 below. In particular for o? = 2 and o? = 3 we have the densities 

pf)(x) = _L (^^2 _ (^2 ^ 27rte + 2iV)e-^*-) , (23) 

pp) (x) = (^48 - (48 + 24tx + Qt^x^ + t"^ x^)e-^^''^ . 

We use now once more (19) but this time with (pk{f) = ^(/) being the j-th intrinsic 

if 



volume of the A;-polytope f,0<j<k. Define (Pki^it)) -. df'"^^ and use [40, Thm 10.3.3] 



to obtain 

^{k,d) ^ r,d-k-i d-k fd-j\fd\ / Kd-i\ '^'^ Kd^a-j 
^ d-j\d-k) \j) \ dKd ) Hj 

as the value for the density of the j-th intrinsic vohme of the (proper) k-skeleton of the 
stationary and isotropic STIT tessellation Y{t). For j = k this yields 

^{k^) _ j{k,d) _ od-k-1 ld\K,d I Kd-1 \ d-k 



g(.,a, ^ ^,K,a, ^ ^d-k-1 r p j td-k^ (24) 

kj Kk \ dKd 
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i.e. Sp"'^ is the mean /c- volume of (proper) /c-dimensional I- faces of Y(t) per unit volume. 
Note, that in particular 

No -N,j-Sy -2 i^a^j^j t 

is the mean number of vertices of Y(t) per unit volume in agreement with (20) above. 
Using the identity 

7V^5.EVolfe(7f ) = 4''") 
combined with (20) and (24) leads to 



for the mean /c- volume of the typical /c-dimensional I- face Especially for A; = 1 and 
A; = d — 1 we have 

EVoMA ) - EVoW(/,_J-^^ j ^. 

In the planar case our formulas specialize to N^^^ — H'^, n[^} — H'^ and Lj — E Voli(/i^^) = 
J and for space dimension d — 3 we reproduce the results from [46], namely A^q^^ = 
(this is the mean number of T- vertices in the sense of [46] per unit volume), N^^j = ^t^, 
N^^j = and Ly — Sy'^^ = Thus, the mean length of the typical I-segment in the 
3-dimensional case equals EVoli(/f^) = f and the mean area of the typical I-face is given 
byEVol2(/f) = ^. 

Theorem 3 may also be used to calculate the mean f-vector {fo''\ f'^\ ■ ■ ■ , f^l^i) of the 
typical /c-dimensional I-face i.e. /j*^'' is the mean number of j-dimensional I-faces of 
the typical k-dimensional I-face ijf^ (thus /q'^^ is the mean number of vertices of a fc-face) 
for < j < k < d, where for k = d we abuse the notation and let ijf"^ stand for the typical 
cell of Y{t). Theorem 3 together with the crucial property (a) of STITS tessellations 
immediately imply that this vector is the same as the mean f-vector of a /c-dimensional 
Poisson polytope. Prom [40, Thm. 10.3.1] we deduce now 



'k' 

In particular, the typical cell of Y{t) has in the mean 2"^ vertices, whereas the zero cell 

has cf. [40, Thm. 10.4.9]. Por the typical cell in the case d = 2 this yields the mean 
f-vector (4, 4) and in the spatial case d = 3 we obtain for the typical 2-dimensional I-facet 
and the typical cell the mean f- vectors (4,4) and (8, 12,6), respectively, which is already 
known from [46]. We have thus proved 
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Corollciry 3 The mean k-volume of the typical k-dimensional I-face is given by 
(,) _ d fdK.Y 1 _ d ( 2^V {^)\' 1 



Voi.(4 ) - _ J - _ 1^ r (I) ) t^ ^^^^ 

and the mean f-vector of if equals {ff' , . . . , /^^\) with ff^ = 2^"^' Q) forO <j <k-l< 
d. 



Beside the typical geometry of the boundary of the typical /c- dimensional I-face, some 
information may also be deduced about the geometry of the relative interior of I^^^^ . Denote 
by N^'^l- the the mean number of (proper) j-dimensional I-faces in the relative interior of 
the typical /c-dimensional I-face if and observe that this mean value is given by 



Nf id-k){t) 



= T7^ = ::;^; 2^"^ o<j<k<d, (26) 



because each such j-dimensional I-face is contained in the relative interior of precisely one 
k-dimensional I-face. For k — d—1, i.e. when regarding I- facets, the formula speciahzes to 

<i-..-^(,')2— ('->-- 

and we obtain the mean value formula 

For example we have in the planar case = 2 in accordance with [32] and for d = 3, 

k = 2 and j = 1 the value A''}^2 i = 4, which fits with the results obtained in [46]. 
Our next goal is a formula for the specific (or mean) j-th intrinsic volume ¥,Vj{lf), < 
j < k, in the sense of [40] of the typical A:- dimensional I-face if . In particular KVk{lf) will 
again be the mean A;-volume, 2EVk-i{lf ) the mean surface area, ^EV,{lf) the mean 
breath and, trivially, 'EVo{lf) = 1 the mean Euler-number of the typical A;-dimensional 
I-face if . Such a formula for the specific intrinsic volumes could easily be obtained by 
dividing their intensities d^'^''^^ by the face intensities Nf . For explanatory reasons we will 
give another proof though, highlighting the temporal construction of STIT tessellations 
and the underlying ideas. 

Corollciry 4 The specific j-th intrinsic volume of the typical k-dimensional I-face of a 
stationary and isotropic STIT tessellation with surface intensity t is given by 

= 7T^ ( : ) I '-"""Wlh^ ) ^7- (28) 
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Proof of Corollary 4 In analogy to what we did in Corollaries 1 and 2, we regard the 
collection of all fc-dimensional I-faces of Y(t) as a marked point process of /c-faces in the 
space of A;-polytopes, where the marks are given by the birth times of the faces and the 
birth time of a k-face is defined to be the birth time of the unique {d — l)-facet which 
induces the birth of the /c-face. Given now the birth time /3 = s of the conditional 
specific j-th intrinsic volume of the typical A;-dimensional I-face ijf^ is 



Kj) \ r (; 



j) \Kk-iJ KjiXksy Kj\jJ [ r(|) 



with Afe as in (3). Here wc have used Crofton's formula for random tessellations [40, Eq. 
(10.28)] and key properties (a) and (b) for stationary and isotropic STIT tessellations. 
Now, Theorem 3 and Corollary 2 say in this context that the the birth time density of ijf' 
is given by 

ds"^-^ ^ 

— < s <t. 

Integration of E,[Vj{ljf^)\P — s] with respect to this birth time density yields the value for 
El^,(/f),i.e. 



d 



fk\ /vn^y 1 

\jJ [ r (f ) j P 




{d-j)Kj\jJ\ r(|) 

and completes the proof. □ 

For j = k, (28) specializes to formula (25), for j = we have the constant value 1, 
which is clear since I^^^ is a convex set with probability 1, and for j = 1 and j = A; — 1 we 
obtain 



kd v/ir(^) i 
d-1 r(|) ? 

A-^/ / 2v^r (^) \ - J 
{d-k + i)Kk-A r(f) i i^- 



k-l 



It is our next goal to generalize some mean value formulas known from earlier papers for 
lower-dimensional stationary STIT tessellations to arbitrary space dimensions. For this 
reason we introduce the notion of k- dimensional J- faces. By such a face for Y{t) we mean 
any /c-face of a (i-dimensional cell, not necessarily maximal. Moreover, a J-face is counted 
with multiplicity given by the number of cells in whose boundary it lies. This way, the 
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collection of A;-dimensional J-faces with such multiplicities coincides with the collection 
of all k- faces of the tessellation cells in Y{t), again with respective multiplicities. This 
convention, as adopted in the existing literature, has the clear and important advantage 
of ensuring that the typical A;- dimensional J-face of Y{t), denoted by jjf^ below, has the 
same distribution as the typical k-face of a corresponding Poisson hyperplane tessellation 
PHT(tA) with the same intensity measure tA, because of the Poisson typical cell property 
(a) of STIT tessellations. Consequently, 

EVjiJlf^) = El^i (Typical k-Face(PHT(tA))), < j < k < d. (29) 

A look at the display (28) in Corollary 4 and a comparison with the second formula on 
page 490 of [40] leads now from (29) to the following comparison relation 

m4'^)-^mjl% o<j<k<d, (30) 

and especially for j = k we have 

EVol,(4'^^) = ^EVol,(jf ). 

Writing now Njfl for the mean number of /c-dimensional J-faces per unit volume we get 
therefore 



E VoU(jf ) . ATg = VoU(4^)) ■ ATg. 



On the other hand, taking into account that each point of the k-skeleton of Y{t) belongs 
to precisely (d — A; + 1) of the k-dimensional J-faces we see that 

E Volfe(4'')) • ATg ^{d-k + l)S'f'''\ 

Putting these together we are led to 

= '-^^^^^t'l (31) 



d 



For d — 2, k — 1 and d — 3, k — lork — 2 these relationships are known from [32] and 
[46]. The value for A^f] can for example be used to obtain a formula for the mean number 
of vertices ^ in the relative interior of the typical /c-dimensional J-face: 

This yields in the planar case | and for d = 3 and k — 2 the value |, both known from 
[32] and [46], respectively. Also observe that the mean f- vector of jjf' is the same as the 
corresponding one for the typical I- face 
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Reiiicirk 1 We wish to remark that the approach used in the proof of Corollary 4 o-nd a 
similar approach used in earlier papers on STIT tessellations, which is essentially based 
on marked point processes, is now fully justified by the considerations of this subsection, 
especially Corollary 2. The potential underlying Slivnyak-type theorem for I- faces and I- 
facets is replaced here by our martingale technique and the fact that (7) is a martingale 
can in some sense be seen as such a Slivnyak-type result for STIT tessellations or more 
generally for iteration infinitely divisible random tessellations. 

3.3 Higher Moments for Typical I-Faces 

This section links mainly to the results from the last subsection and it is our goal to establish 
formulas for the second and third moment of the fc-volume of the typical fc-dimensional I- 
face of stationary and isotropic STIT tessellations Y(t) with surface intensity t. Moreover, 
we will establish a formula for the second moment of the {k — l)-dimensional surface area 
of the boundary of I^''^ . For the special case k — 1, formulas for the second or higher 
moments of the length of the typical I-segment can be obtained directly from the explicit 
density (22). For l<A;<d — Iwe calculate 

EVol^(/;'^^) = / / Xise-^''''^-dsdx= / dXi^ x"" e'^'''' dxds 
Jo Jo t Jo t Jq 

^ d.n\ d-n\ ( V^r{^) y 1 

t'^X'l Jo d-ny r{^) J t""' ^ ' 

For k > 1, i.e. for higher dimensional I-faces, we start with the following general state- 
ment, which directly follows from Theorem 3 and the fact that TypicalCell(PHTfe(Afes)) 
has the same distribution as the typical /c-face Typical k-Face(PHT(s)) of a stationary 
and isotropic Poisson hyperplane tessellation with surface intensity s in MJ^, where by 
TypicalCell(PHTfc(Ajfcs)) we mean the typical cell of a stationary and isotropic Poisson 
hyperplane tessellation in M'^ with surface intensity A^s and Xk given by (3). 

Corollary 5 The n-th moment of the j-th intrinsic volume of the typical k-dimensional 
I-face ijf^ of a stationary and isotropic STIT tessellation Y{t) in M'^ with surface intensity 
t satisfies 

ET//(/f^) = f k-Face(PHT(s)))ds 
Jo i 




EV^/(TypicalCeU(PHTfe(Afes)))cis (34) 



with 0<j<k<d— 1. This value is finite if and only if d — jn > 0. 

At first, we apply (34) for j = k and n = 2. It is well known that in this case 

EVol2(Typ,calCell(PHT.(A..))) . |(^)^'^ . | (|!l)''-L, 
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see [22, p. 179] and we obtain by integration 



^^^^'^^^'^ ^-^^2^1^ r(|) J ^' 

whenever d — 2k > and +oo in the other cases. For d >S and A; = 1 our formula yields 
the second moment of the length of the typical I-segment, see (33), which speciahzes for 

= 3 to the value E Voli(/f ^) = f|. 
In view of the results from [22, Sec. 6.3], a similar formula is also available for the third 
moment of Volfc(/j^'^'*), i.e. for the case j — k and n = 3. For c? — 3A; > we have 



3fc 



EVol^(/f ) = :r^2-7r-^ 



d-3k V ^(^) / r(^5^) 

and E Vol^(/^'^'*) = +oo, whenever d — 3k < 0. In particular for d > 4 wc obtain once more 
the third moment of the length of the typical I-segment (33) with n — 3. By the same 
method we also get for d > 3 and k > 2 a formula for the second moment of the surface 
area of the boundary bd(/^'^^) of ijf^ (this corresponds up to a factor 2 to the case j — k — 1 
and n = 2 in Corollary 5): 

EVolti(bd(/f )) = 




if d + 2 — 2A; > 0, see again [22]. Especially for d = 3 and k — 2 we obtain the second 
moment of the perimeter of the typical /-facet 

EVol?(bd(/f)) = ^('l+"'^ ' 



4 V 4 y ^2 

For the particular case d = 3 and k = 2 we also know that the second moment of the 
number of vertices of equals 

Explicit formulas for other or higher moments of the intrinsic volumes of ijf^ are not 
possible, since the corresponding values for Poisson hyperplane tessellations are currently 
not known up to our best knowledge. 
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3.4 Typical Geometry beyond the Isotropic Regime 



It is our aim to explore in this section some aspects of the typical geometry of STIT 
tessellations and more generally iteration infinitely divisible random tessellations beyond 
the isotropic regime considered in the last paragraphs. We first put ourself in a very general 
position by taking A to be an arbitrary non-atomic and locally finite measure on the space 
H of hypcrplanes, which is absolutely continuous with respect to the motion-invariant 
measure Aiso- In this case the resulting tessellation Y{t), regarded as a particle process of 
convex cells in the sense of [40, Chap. 4], admits a translation-regular intensity measure 
in the sense of [40, Chap. 11] and the theory presented there applies. Furthermore, 
the intensity function Njf'j(z), 2; £ M*^ of /c- dimensional I-faces, the intensity function 

d!^l^\z), z ^ oi the j-th intrinsic volume of the /c-skeleton and the lower-dimensional 
directional distributions in the sense of [40, Chap. 11] of the iteration infinitely divisible 
random tessellation Y{tK) are well defined. Combining the results from [40, Chap. 4 and 
Chap. 11] with Theorem 2 from above yields at first 

Corollary 6 The directional distribution of the random process of k-dimensional I-faces of 
an iteration infinitely divisible random tessellation with translation-regular intensity mea- 
sure A coincides with the directional distribution of the {d — k)-th intersection process of a 
Poisson hyperplane tessellation with the same intensity measure. 

In the translation-invariant case and only for d = 2 this was already shown in [24] by a 
quite different argument. 

Denote by 7pHT(sA) (^) intensity function of fc-faces of a Poisson hyperplane tessellation 
with translation-regular intensity measure sA. Then we obtain by applying Theorem 2 
above and [40, Thm. 11.4.2] the local mean value formula 

N'Siz) = (rf - A:)2"-'-' £ lj^iT(.AM>ds = (<i - k)2-->'-' ^' lj'°l^,.^^(z)ds 

for Lebesgue almost all z E M!^. Analogously, for the intensity function d^p{z) of the j-th 
intrinsic volume of the /c-skeleton of Y{tK) we obtain the local mean value relation 

for almost all z E wrt. Lebesgue measure. In the translation-invariant set-up, the 
considered intensity functions are constant and the mean value formulas hold without 
reference to the points z G W^. Beside the two mean value relations from above, general 
local mean value relation may be obtained by generalizing (19) to the non-stationary case 
by a suitable concept of localization, see [40, pp. 523-524]. However, the absence of scaling 
relations in the non-stationary case makes them less explicit. 
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We pass now to the stationary but not necessarily isotropic regime and assume for the 
remaining part of this subsection that the measure A is translation-invariant rather than 
just translation-regular. In this case, the result of Theorem 3 remains also true, with 
PHT(s) replaced by PHT(sA) there. To see it, observe that the mean value formula 

'■^ ~ d 'PHT(tA) 

is in force also in the stationary but non-isotropic case, in full analogy to (21). Moreover, by 
translation invariance there exists some constant c, not depending of such that 7pHT(tA) ~ 
cf^. Thus 

1 7PHT(.A) ^ 1 cs'^ ^ d s'^-' 

s jyi^l ~ s d^2'^->'-^ct'^ ~ {d- k)2'^-^-^ 

and we have again in analogy to (21) 

(fc) 



and hence 



l'^Q'r"'"'ls (36) 



Y(t) ■ '■"'^'^ 



also in the non-isotropic case, where by Q™^*^**^'' we understand the distribution of the 
typical /c-dimensional I- face of the stationary Poisson hypcrplane tessellation with intensity 
measure sA and where Y{t) — Y{tA) is the STIT tessellation with intensity measure tA. 
Also the time-marked analogue of (35) holds true in the non-isotropic setting: 



* ds'^ ^ ^PHT(sA) ^ r 
Ik 







ds. (36) 



Note that the statement about the mean f-vector of the typical A;- dimensional 1-face in 
Corollary 3 remains valid as well as the mean value formulas (26), (27), (30), (31) and 
(32). 

The equality (36) from above allows us to determine the conditional distribution of the 
typical A;-dimensional 1-face of Y{t) = Y{tA) given its birth time < s < t. We readily 
sec that this conditional distribution equals the distribution of the typical fc-face of a 
Poisson hyperplanc tessellation with the same parameter measure tA. Especially, the 
conditional length distribution of the typical I-segment given its birth time s and direction 
I is an exponential distribution with parameter sA([e(/)]) where e(/) is a unit vector on 
I, and thus the distribution of the length of the typical I-segment itself is a mixture of 
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suitable exponential distributions. In particular, our considerations answer an extended 
version of the question formulated in the Outlook section of [24] . Moreover, the directional 
distribution 71^^^ of the random process of /c-dimensional I-faces, 1 < /c < ci — 2, is in view 
of Corollary 6 and (4.62,4.63) in [40] given by 

where by p{d-k) we denote the {d — k)-th projection generating measure of the associated 
zonoid n [40, (4.46)] - an auxiliary convex body, which uniquely characterizes the law of 
the underlying STIT tessellation - with generating hyperplane measure tA, as defined by 
[40, Eq. (14.36)], whereas (•)-*- stands for the orthogonal complement mapping. 
Beside the distributional result for I-segments, mean values for the typical /c-dimensional I- 
face of a stationary but non-isotropic STIT tessellation arc also within reach of our methods 
and can be obtained by the same scheme as demonstrated above. In place of the concrete 
values involving Gamma-functions, in the non-isotropic case we have at our disposal the 
geometric parameters of the associated zonoid. As above, let 11 be the zonoid in the sense 
of [40, Chap. 4.6] with translation-invariant generating hyperplane measure tA. Then the 
j-th intrinsic volume of the typical /c-dimensional I-face is given by 

' ' ~ Jo 0voi4fn)''- Q voun)io (f)^ 

Q d-jVohiU)- ^ > 

Here we have used (35), [40, Thm. 10.3.3] and the homogeneity of the intrinsic volumes. 
Note that (37) reduces to (28) in the isotropic case, i.e. when A = Aj^o or equivalently 
when the zonoid 11 is a d-dimensional ball with radius proportional to t. 

4 Second-Order Properties 

In this section we study second-order characteristics of iteration infinitely divisible ran- 
dom tessellations and stationary STIT tessellations. This is done from several different 
perspectives, both with the use of our martingale techniques and of recent very elegant 
second-order theory developed by Weiss, Ohser and Nagel [47] for the planar case. 

4.1 Martingale Tools 

The general martingale statements of Section 2 admit a convenient speciahzation to deal 
with the second-order characteristics of iteration infinitely divisible or stationary STIT 
tessellations as well. Taking to be a general bounded measurable functional of {d — 1)- 
dimensional facets, regarded as usual as closed subsets of W, we put G{Y,t) := {T,^{Y) — 
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EE^{Y{t, W))y so that G{Y{t, W),t) ^ ^iYit, W)) with E^(y(i, W)) := E^(y(i, W)) - 
EE<^(y(i, W)), and we use (10) to check that 

§^G{Y{t, W),t)^ -2[E^{Y{t, W)) - EE^{Y{t, W))]EA^{Y{t, W)) (38) 



where 



A^{Y):^ [ V (j){Facet{x;H\Y))A{dH). (39) 

Put together (6), (8) and (38) and use locahzation as in the discussion preceding (10) with 

Gn, N ^ oo, chosen so that (Gjv(-, •) A A^) V -A^ = {G{Y,t) A N) y -N, that \Gn{-, •)! < 
+ 1 and that Giv(-,t) be twice continuously differentiable in t, and with the locahzing 
stopping times = mit>o{\G{Y{t,W),t)\ V \§^G{Y{t,W),t)\ V \^G{Y{t,W),t)\) > N. 
Proceeding as there, we readily conclude that 



S2(F(t,iy))- / / (j)'^{Fiicet{x;H\Y{s,W)))A{dH)ds+ 

hi 0(Facet(x;//|y(s,iy)))[E^(y(s,iy))-EE^(y(s,iy))]A(di/)- 
[E^(y (s, W)) - EE^(y (s, W))]EA^{Y{s, W))]ds = 

tl{Y{t,W))- f A^.{Y{s,W))ds-2 f A^{Y{s,W))E^{Y{s,W))ds (40) 
Jo Jo 

is a ^^t-martingale as well, with A^{Y{s,W)) := A^{Y{s,W)) - EA^{Y{s,W)). Note for 
future reference that taking another bounded measurable facet functional ip, applying the 
above for + and 4> — ip and subtracting yields one further Q^^-martingale 

s^(y(t,iy))s^(y(t,iy))- / A^4Y{s,w))ds- 

Jo 

f {A^{Y{s,W))f:4Y{s,W)) + A4Y{s,W))i:^{Y{s,W)))ds. (41) 

For general this cannot be simplified any further. However, in our further considerations 
we shall focus our attention on translation-invariant face functionals of the form 

0(/):=Vol,_i(/)C(n(/)) (42) 
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with n(/) standing for the unit normal to / and C for a bounded measurable function on 
<S(i_i. Then, using (42) we see that 

A^= [ Vold-i{HnW)C{n{H))A{dH)^ const 

J[W] 

and so ^4,^ = and thus, by (10) and (40), 

E^{Y{t,W)) and mY{t,W))- [ A^2{Y{s,W))ds (43) 

are both $5t-martingales. In particular, see [16, Thm. 4.2], the martingale Tj^{Y{t, W)) has 
its predictable quadratic variation process {T,^{Y{-, W))) absolutely continuous and given 

by 

{j:^(Y(;W)))t^ f A^.{Y{s,W))ds. (44) 
Jo 

These observations arc going to be crucial for our second-order analysis of iteration in- 
finitely divisible and stationary STIT tessellations, given in the next Subsection 4.2. 

4.2 Variance of the Total Surface Area in a Window W 

It is the main purpose of this subsection to calculate the variance of the total surface 
area of the iteration infinitely divisible random tessellation Y{t, W), for a compact convex 
window C R'^ with the property that Vold(l^) > 0. We let A be an arbitrary diffuse 
and locally finite measure on T-L. Recall now (39) and note that it implies 



A^.{Y{t,W))= [ Yl <t>\mdH) 
/eCeiis(ynH) 

/ C^(n(i/)) / / y are in the same cell of F n H]dxdyA{dH). (45) 

Jiw] JwnH JwnH 



l[W] 

Thus, using (43) and taking expectations of both sides yields immediately 

VarS^(F(t, W)) = 

C^(n(i/)) / / P(x, y are in the same cell of Y{s, W) n H)dxdyK{dH)ds. 
J [w\ JHnw JHnw (-45-) 

Taking into account that 

F{x, y are in the same cell of Y{s, W) H H) = exp{—sA{[xy])), 
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which follows from key property (a) and using (46) we end up with 

Var(E^(y(i,iy)))= / / C\n{H)) [ [ exp{-sA{[xy]))dxdyAidH)ds 
Jo Jiw] JwnH JwnH 



'0 J[W] JwnH JwnH 

= / (Hfimf f i^i^||z«ldxd,A(.7/). (47) 

Jiw] JwnH JwnH ^^[[^y\) 

For the stationary and isotropic case A = Aj^o we want to evaluate this integral further in 
the special case = Vol^-i, i.e. when ( = 1. First, we use the affine Blaschke-Petkantschin 
formula [40, Thm. 7.2.7] with g = 1 to obtain for any non-negative measurable function 
h : (R"')2 ^ E 

h{x,y)dxdy ^ ^ J J J h{x,y) \\x - y\f~^ eL{dx)eL{dy)dL, 

where C is the space of lines in M.'^ with invariant measure dL, i.e. the affine 1-dimensional 
Grassmannian in R'^ and £l is the the Lebesgue measure on L with normalization as 
specified in [40, Thm. 13.2.12]. Taking now 

h{x,y) = l[x e W]l[y e W] \\x-yfg{x,y) 

for some k > —d and another non-negative measurable function g : {W^Y we obtain 

X -yt g{x,y)dxdy j [ [ \\x - y\\'^'^^'' g{x,y)iL{dx)diL{dy)dL. 

wJw ^ JjC JwnL JwnL ^4g^ 

For k = —1 this yields 

'^\dxdy = II I \\x — yW^ ^ g{x,y)iL{dx)iL{dy)dL. 

'c JwnL JwnL ^49^ 



w Jw \\x — y\\ 2 



We replace now in (48) for k = 0, W hy W (1 H for some fixed hyperplane H and d by 
d — 1 and get 

/ / g{x,y)dxdy 
JwnH JwnH 

= (llll)^ II I _ yf-^ ^(^, y)e4dx)e4dy)dL^, 

^ J c" JwnHnL JwnHnL 

where by we mean the 1-dimensional affine Grassmannian restricted to H with invariant 
measure dL^ . Averaging the last expression over all hyperplanes H and using the fact that 



29 




H JwnH JwnH 



Aiso{dH) (8) dL" = dL, see [40, Thm. 13.2.12], yields 

g{x,y)dxdyAiso{dH) 

= - 1)^^-1 f II I \\x-y\\''-^g{x,y)^L{dx)^L{dy)dL^Kso{dH) 

^ Jh J ch JwnHnL JwnHnL 
^ {d- l)^,_i I I I ||^_^||<^-2^(^^^)^^(^^)^^(^^)^^_ (50) 

^ Jc JwnL JvvnL 
By comparing (49) and (50) we finally conclude the non-trivial identity 

/ / / g{x,y)dxdyAUdH)^^^-^^ f f-f^My (51) 
Jiw] JwnH JwnH ^>^d JwJw\\''^~y\\ 

for any non-negative measurable function g : {W^Y — >■ R. Using Equation (14) in [31] 
together with the the mean projection formula [40], Thm. 6.2.2 for q — j = d — 1 (or 
equivalently using the sectional property (b) in Subsection 1.2) we get from (51) with 

1 - exp(-tA,,o([xi/])) 1 - e-^*l'"-^ll 



AisoiXxy]) \\x - y 



dKa 

the following formula for Var(Evoid_i(l"(i, W))) - Var(Vold_i(r(t, W))): 



f f f d— Iffl — e'''^d, ^" 

/ / / 9{x,y)dxdyKiso{dH) ^ —— I I dxdy 

J[w]JHnwJHnw ^ Jw Jw \\x — y\\ 



2k 



d{d - l)Kd 



lwiry~^ ( 1 - e"^^*M dr 



2 

by using li- dimensional spherical coordinates. Here 



J Sd-i 



is the isotropized set-covariance function of the window W . Summarizing, we arrive at 

Theorem 4 For the stationary and isotropic STIT tessellation Y{t) with surface intensity 
t > we have 

-l\\:r-y\\ 



Var(Volrf_i(F(t,W^))) = ^— / / ^ dxdy ^ (52) 

^ JwJw \\x-y\\ 

T.M.-(l-e-^'')*. (53) 
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We specialize this now by taking W to be the ball in R'' with radius it! > 0, i.e. W — Bj^. 
In this case, the isotropized set-covariance function (r) is given by 

7.. (.) ^ 2«^«._. /V - u^)^,u ^ 2..«^ (i - (i, i^; I ^)) 

(with 2-^1 being the Gauss hypergeometric function) for < r < 2R and otherwise, see 
[11, Chap. 4.8.4], and we have 

Var(Vold_i(F(t, 5^))) = d{d - l)KdKd-iR'^ / ( 1 - e"^^*" ) / ~ u^)^dur'^-^dr 

rfi-e-^'-l i^;|;^lV'-'*. 



2 2R Kd V2' 2 ' 2' 4i?2 
where the constant before the integral may also be written as 

d{d - l)nand-,R' = (d - 1)^-^^^^R' = ^-^^R' 

with Ai from (3). In the even more special case d — ?> the isotropized set-covariance 
function 7^3 (r) takes the form 



f^'(l-|i + Tt) :0<r<2i? 

: r > 2i? 



and the variance integral can be evaluated in a closed form: 

Var(Vol2(F(t, = (i2i?2(i2 _ stR + Si^i?^) + 24(1 + ii^je"*^ - 24) . 

(54) 

The same closed form cannot be obtained for d = 2, since 7^2 (r) has a more complicated 
structure, i.e. 



7^2 (r) = 2R^ arccos (^) - ^V4i?2 _ 



for r between and 2i? and 7^2 (r) = for r > 2i?. Unfortunately, the resulting integral 
can in this case not further be simplified. 

Another important task in this context is to determine the large R fixed t asymptotics of the 
variance Var(Vold_i(F(t), Wr)) for the family of growing windows Wr = R-W, i? — )■ 00. 
For d — 2 we claim that 

Yar{Yoh{Y{t,WR)) - tt Vol2(IF)i?^ logi?, (55) 

where ~ stands for the asymptotic equivalence of functions, i.e. f{R) ~ g{R) iff f{R)/g{R) ~^ 
1 as i? — >■ 00. Indeed, this can be estabhshed by using (53), the relation 7^^^ ~ Vol2(iyR) = 
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R'^YohiW) valid uniformly for argument r = 0{R/\ogR), the observation that — > 

for r = Q.{R\ogR), together with the fact that — e~'^^)y ~ log-R, c > 0, as soon 

as log L(i?) ~ logi?, and the key property (c) of STIT tessellations: 



VaT{yoh{Y{t,WR))) = r2Var(Voli(r(l,mi?))) = ^^ fw^'^-e-^- 



dr 
r 

nr^ YohiWtR) log(tR) = TT Vol2(M^)i?^(log R + log t) 
7rR^Yo\2(W)\ogR. 



Another way to see it is to combine (63) below with the results from Subsection 5.2. Thus, 
as we see, in the planar case the studied asymptotics only depends on the area of W. Things 
get more comphcated for d > 2 though. To see it, use (52) and the scahng property (c) of 
STIT tessellations to obtain 



yaT{Vo\a-i{Y{t,WR))) = i?2('^-^)Var(Volrf_i(r(i?t,iy)) 



2 JwJw \\x-y\\ 



Consequently, we get for d > 2 and — > 00 



Var(Vol,_i(r(i, Wn))) - R'('^-^)i-lE,{W) (56) 



where £'2(1^) is the 2-energy ofW, see [23, Chap. 8] given by 



E2{W)^ I I \\x-y\\-^dxdy. (57) 

Observe that this docs not extend for the separately treated case d = 2 because there the 
integral in (57) diverges. It is easily seen that E2{-) enjoys a superadditivity property 

E2{Wi U W2) > E2{Wi) + E2{W2), W^i n = (58) 

which stands in contrast to (55) where the asymptotic expression is linear in Vol2(W^). 
We will now derive an integral geometric expression for this energy functional. Taking 
g(x, y) = 1 and A; = — 2 in (48) yields the remarkable identity 

MW) = / / \\x-y\\-'dxdy^^ [ [ [ \\x - y\f-' eUdx)hidy)dL 

- a - m - 2) /.^"'-'^ ^ ^'"'^^ ^ - m - 2) ^"''^'' 

with Id-iiyV) being the {d — l)-st chord power integral of W in the sense of [40, p. 363]. To 
display the dependency of the asymptotic variance Var(Vold-i(F(t, Wr})) on the geometry 
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of W we could take instead of the {d — l)-st chord power integral Id-i{W) of W also the 
functional 

vl,(wnE)dE, 



iA(d,d~2) 

where A{d, d — 2) is the set of affine {d — 2)-planes in W'- with invariant measure dE, since 

h-AW) = ^A^^ [ vUw n E)dE 

^l^d-l JA{d,d-2) 

according to [40, Eq. (8.57)]. Thus, comparison with (59) yields the alternative represen- 
tation 

E,(W) = , , "^^l I vUw n E)dE. 

The choice of one of these functionals is more or less a matter of taste, but we will take 
Id-iiW) here, because it will allow us in Section 6 a better comparison with other tessel- 
lations models. Hence, combining (59) with (56) from above, we arrive for d > 3, i? — > oo 
at 

Var(Vold_i(r(t, Wr))) ~ -^I^_^{Wr) = -^B?i'^-^) I^_^{W). (60) 

a — z a — z 

In general, la-iiW) cannot further be evaluated. But for W = Bf we have by applying 
[40], Theorem 8.6.6 (with a corrected constant) 

I,_,{Bi)^d2'-^^:^!^ 
i^d-i 

and, thus, the 2-energy of the ti-dimensional unit ball Bf equals 



^'^^'^ {d-l){d-2) Id-i (d-l)(d-2)^[2 

In the particular case ci = 3 we obtain the value E2{B\) = 47r^, which agrees with the 
explicit variance formula (54). Another case, where the chord power integral can be eval- 
uated is the practically relevant case of the cube in with edge length a > 0. Here 
we have 

57r 13 1 
I^{Cl) = — - — - 2\/27r + 27rln 2 - - In 2 + 4^2 arctan \/2 - IQF ^ 3.7557 
3 9 3 

with 

arctan s/x'^-2 



dx, 

for a = 1 and in general hiC^) — 0^/2 (C*?) which can easily be obtained from [15, Eq. 
(15)] and the fact that Id-i is homogeneous of degree 2d — 2. 
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Remcirk 2 Note that in the 2-dimensional case the asymptotic variance of the total edge 
length in Wr may also be written as 

Var(Voli(r(t, W^i?))) ~ Ii{W)RHogR, 

yielding some kind of consistency with the higher dimensional cases represented by (60). 
But the argument using 2-energies given for d > 3 above cannot be applied to prove this 
formula. 

Remcirk 3 In our situation, [40, Thm. 8.6.5] for 

f(r)^f(\\x-y\\)^^^^g{x,y) 

11*^ y\\ 

can be applied, to deduce that Var(Vold_i(F(t, W))) is maximal among all compact convex 
bodies W with positive volume exactly for d-dimensional balls, i.e. when W — Bj^. 

Remark 4 In particular, Theorem 4 establishes weak long range dependence (see [39]) 
present in stationary and isotropic STIT tessellations Y{t), since 

Var(Volrf_i(y(t,PyR))) _ \c2\ogR^oo :d^2 
YoUWr) ^ \cdE2{W) Vold{W)-^R'^-^ ^ oo : d > 3, 

as -R —)■ oo with being space dimension dependent constants. Heuristically this can be 
explained by the geometry of the I-facets ofY{t). It was shown above that the [d—\)-volume 
of the typical I-facet has first but no second moment (see [25] in the planar case or [45] 
for a much more detailed discussion). Thus, the I-facets ofY{t) can be extremely large, 
but there are no 'full' planes in the tessellations. This is reflected in the planar case by 
the log-term in the asymptotic expression for the variance of the total edge length, whereas 
in higher dimensions by the non-additivity of 2-energies determining the asymptotics. We 
will compare this with other tessellation models at the end of this paper. 

Remark 5 In the planar case, for W = B\ the unit ball andt — 1, (55) was independently 
established by Lothar Heinrich (personal communication) using a quite different method 
which allows also the derivation of (55) in its general form in the planar case d = 2. The 
main idea is to dilate Y(t, Bj^) - regarded as a random closed set in the plane - by a small 
ball with radius e and to use the capacity functional ofY{t,Bjf) (see [31]) together with a 
generalized Steiner formula and a limiting argument. Unfortunately the asymptotic formula 
for the variance is not sufficient to establish a central limit theorem as will be shown below 
in Section 5. 
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4.3 Pair-Correlation Function and an Alternative Approach in 
the Planar Case 

The second moment measure iJ,^^^ (■><■) of the random surface area measure of the stationary 
and isotropic STIT tessellation Y{t) C can be defined by 

f^f{A, X A2) := E[yoU-i{Y{t,A,))yoU-i{Y{t,A2))], 

where Ai and A2 are bounded Borel sets in R''. By translation-invariance of the surface 
area measure of Y{t) it is sufficient to regard the the reduced second moment measure fCd 
of the random surface area measure of Y (t) , see [43] , and since we are also in the isotropic 
case it is even enough to consider the reduced second moment function Kd{r) given by 

Kd{r):^K:d{B^), r > 0. 

Then quantity t ■ Kd{r) can be interpreted as the mean surface area of faces of Y{t) within 
a ball with radius r centered at a typical point of Y{t) (when the tessellation is regarded 
under a suitable Palm distribution). In the case where Kd{r) is differentiable in r we can 
consider the pair- correlation function gd{r) given by 

1 dKdjr) 
''^'^ ^ d.dr^-^ dr ' 

It describes the expected length density of Y{t) at a given distance r from a typical point 
oiY{t). 

Prom the variance formula in Theorem 4 the following can be easily deduced: 

Corollary 7 The pair- correlation function gd{r) of the random surface area measure of 
the stationary and isotropic random STIT tessellation Y{t) is given by 



9d{r) = 1 + ^ I 1 - e 



2k 



dKj 



Especially for d — 2, gd{r) becomes 

92{r) = 



2t^r 



(l-e-^*), 



which was independently obtained by V. Weiss, J. Ohser and W. Nagel and will be pre- 
sented in the forthcoming paper [47]. We want to briefly sketch an alternative approach 
to obtain the variance of the total edge length of Y{t,W), d = 2, by using their recent 
results. The main result of [47] says that the second moment measure of the random length 
measure induced by Y (t) coincides with the second moment measure of the length measure 
induced by a Boolean segment process whose primary grain has length density pi{x) given 
by (23), which is the density of the length distribution of the typical I-scgment of Y{t) in 
the planar case. Por a Boolean segment process with length distribution function L and 
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mean I of the typical segment we know that its reduced second moment function K2{r) is 
given by 

K2(r) = Trr^ + i ( /" x^L{dx) + I {2xr - r^)L(dx)) , 

LaI \Jo Jr J 

see [42]. Here is the mean segment length per unit area. In our special case, with 
1 — J and La — t and in view of the fact that length density (23) is a mixture of densities 
of exponential distributions, where the mixing density equals ||, by applying repeated 
integration by parts we are led to 

K2{r) = Trr^ + ^ + In 0tr^ + Ei 0ir^ ^ , (61) 

with 7 being the Euler-Mascheroni constant and Ei the exponential integral 

/•oo ^-xy 

Ei(x) = / dy. 

Ji y 

Using now the definition of the pair-correlation function g2{r) and the fact that the deriva- 
tive of Ei(a;) equals —e~^/x we get the following formula for the pair-correlation function 
g2{T) of the random length measure of Y{t): 



and this is the same as we obtained above in Corollary 7 for d = 2. Using the pair- 
correlation function calculated in this way together with the formula 

POO 

Var(Voli(r(t, W))) = 2TTt^ / ^wir)g2ir)rdr - t'^YolKW) 

Jo 

from [43, p. 233] we can also obtain (53) for d = 2. It should be emphasized though that 
our original approach developed in Subsection 4.2 above yields information also on higher 
dimensional cases. For example we have for the spatial case d — 3 



9sir) = 1 + ^ (l - e-^*^) , 



Ks{r) = l^r^ + |L (3tr - 6 + 6e-i*'-) . 



The general expression for Kd{r) is rather complicated and involves again special functions 
and is for this reason omitted. 



5 Central Limit Problem 

Having studied the first- and second-order properties of iteration infinitely divisible and 
stationary STIT tessellations, we now pass to the central limit problem in the stationary 
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(thus STIT) regime. This problem will be considered in two closely related settings, in- 
terestingly leading to results of very different qualitative natures. First, in Subsection 5.1 
we shall focus our interest on residual length/surface increment processes arising respec- 
tively as cumulative length or surface area of I-facets born after a certain fixed time in the 
MNW-construction. In this set-up we shall establish a functional central limit theorem 
with the limit process identified as a suitable time-change of the standard Wiener process. 
Next, we shall pass to the full length/surface process, taking into account also the I-facets 
born at the very initial big bang stages of the MNW-construction, as descriptively termed 
in [26]. It turns out that whereas in dimension 2 the functional central limit theorem and 
Brownian convergence is preserved as shown in Subsection 5.2, this is no more the case 
for dimensions 3 and higher, where non-Gaussian limits arise as argued in Subsection 5.3. 
This apparently surprising phenomenon is in fact due to the influence of the big bang phase 
itself, which is negligible in two dimensions but turns out crucial in higher dimensions. 
To proceed with our discussion, we put Wr := RW for i? > and with W standing for 
a compact convex set of non-empty interior, to remain fixed throughout this section. We 
also let A be some translation-invariant measure on as in Section 1.2. Finally, we only 
shall consider of the form (42) in this section. 



5.1 Brownian Convergence for Length/Surface Increment Pro- 
cesses 

For fixed sq > we consider the centered surface increment process 



^[f:^{Y{t,Wn))-f^^{Y{so,Wn))] 



t&[so,l] 



corresponding to the increments of S^(y(t, Wr)) relative to some initial time moment Sq 
in the MNW-construction. The main theorem of this subsection is 



Theorem 5 For each Sq > the centered surface increment process 



te[so,i] 



converges in law, as i? —)■ oo, on the space V[so, 1] of right continuous functions with left- 
hand limits (cddldg) on [sq, 1] endowed with the usual Skorokhod topology [2, Chap. 3, Sec. 
14], to time-changed Wiener process 

where >V(.) is the standard Wiener process and Vi^(0, A) is given by (63) or alternatively 
(72) below. In particular, 
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converges in law to J\f{0, Vw{(f), A) J^^ "^ds), a normal distribution with mean and vari- 
ance A) f^^ s^~'^ds. 

Note that this statement cannot be extended to Sq I 0, as would be of interest as poten- 
tially yielding Gaussian limit for T,^{Y{1, Wr)). The problem is that the variance integral 
A) J*^ s^~'^ds diverges at 0. As already signalled above, we will see below that this 
difficulty can be overcome for d = 2 but not for d > 2 where we will obtain non-Gaussian 
limits for suitably normahzed T,^{Y{1, Wr)). 



Proof of Theorem 5 Note first that 
1 



U,.(y(l,iy«)) = i / ^C'inm Yl yolU{f)A{dH) 



I R^rC'lnl^)) E Vol^_i(/)A(d//). (62) 

We claim that upon letting it! — > oo this converges in probabihty to 

. .X f o,^, , ,EVol^ i(TypicalCell(i/ nr(l))) . , , 



where is the orthogonal complement of w e Sa-i and TZ is the directional distribution 

of the stationary STIT tessellations Y{t) as given in (2). To see it, recall that RH fl Y{1) 
is a STIT tessellation in RH for each R > and H E Ti. Thus, applying [40, (4.6) and 
Thm 4.1.3] and (10.4) ibidem to this tessellation, we have 

lim -i-E V Vol^ ,(/) 

/e Ceiis(iiif ny ( 1 , Wji) ) 



VI <^r. ,,,, E\'oli,(Tyi)icalC(41(/f n y(l))) 

= Voi._, n w) EVoi,_,(iypicaiCeii(iiny(i))) • 

Next, we observe that RH n Y(\,Wr) = R-h {H nW)n Y{1) where -h is the scalar 
multiplication relative in if, that is to say H 3 R -h x = ph{0) + R{x — Ph{^)), x E H 
with ph standing for the orthogonal projection on H. Thus, using the strong mixing and 
tail triviality theory for STIT tessellations recently developed by Lachieze-Rey [21, Thm 
2], noting that tail trivial stationary processes are ergodic [10, Prop. 14.9] and then 
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using the standard multidimensional ergodic theorem, see e.g. Cor. 14.A5 ibidem, to 
fli^E/6Ceiis(ii.«(/fniy)ny(i)) Vold_i(/), we get from (64) that 



lim V Vol^ i(/) 

/eCeiis(R//ny(i,H/ij)) 



VI ' VolL,(TypicalCell(ff n ^(1))) 

= Vol,_,(i/ n ^O jEVol._,(TypicalCell(ifnr(l))) ^^'^ 

in probability. Putting this together with (62) and integrating over [1^] yields 

lim ^A^2 (F(l, Wr)) = Vw{(p, A) in probability (66) 

-R-^oo R"- 

as required. Note now that by scaling properties of Y{s, Wr) (see key property (c)) and 
0^ for s > we have 

^,A,.{Y{s,Wr)) ^ ^j^,A,4Y{1,Wsr)) ^ ^j^^A,4Yil,WsR)). 

^ > (67) 

Thus, combining (66) with the scahng relation (67) we get 

lim —lAxiiYis^WR)) — — T— i-lV(0, A) in probability uniformly in s e [sq, 1]. 

(68) 

This crucial statement puts us now in context of the general martingale limit theory. 
Indeed, using (43) we see that = -^[E^(r(l, Wr)) -Jl^iY {sq, Wr))] is a martingale 
with absolutely continuous predictable quadratic variation process 

{Sf^)t = f Y^^Ayi^^ WR))ds, (69) 

see [16, Thm. 4.2]. In these terms, (68) yields for each t 

hm {S^^f), = /* -^Vw(<l>, A) in probability. (70) 



We want now to apply the martingale functional limit theorem, see e.g. [48, Thm 2.1]. 
The condition (ii.6) there is just (70) whereas condition (ii.4) there is trivially verified 
because the predictable quadratic variation (5^'!^) has no jumps by (69). It remains to 
check the condition (ii.5) ibidem, which is that the second moment of the maximum jump 
^{S^^J^; 1) of the process {S^^'^)se[so,i] goes to as — > oo. To this end, note first that, 
with probability one, i7(iS^'^; 1) is bounded above by a constant multiple of i?""^/^ times 
the {d — l)-th power of the diameter of the largest cell of Y{so,Wr). Since the typical 
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cell of Y{sq) is Poisson with intensity measure sqA by property (a) of STIT tessellations, 
we conclude by standard properties that the expected number of cells in Y{sq,Wr) with 
diameters exceeding D is of the order 0{R'^ exp(— D)) since A has been assumed to have its 
support spanning the whole of M'^, see Subsection 1.2. Summarizing, setting u = D'^~^R~^^^ 
we are led to 

F{J{S^^^; 1) > u) = 0(i?'^exp(-i?'^/(2d-2)^i/(d-i)^)^)_ ^^-^^ 
Clearly, (71) is much more than enough to guarantee that 

hm Ej2(5iJf;l) = 

which gives the required condition (ii.5) of Theorem 2.1 in [48]. Upon a trivial time change, 
this theorem yields now the functional convergence in law as stated in our Theorem 5. This 
completes the proof. □ 

An alternative formula can be provided for the factor VV(0, A). Denote by 11 the as- 
sociated zonoid of STIT tessellation with generating hyperplane measure A, by 11° its dual 
body and by TZ the directional distribution of the STIT tessellation. Then we have 

Proposition 1 It holds 

where Yl\u-^ stands for the orthogonal projection o/II onto the hyperplane and where the 
polar body {Jl\u-^)° is considered relative to w^. In the isotropic case, i.e. when TZ — v^-i, 
this reduces to 

Vw{ct>Aiso) = Vol,(W^)2<^-V'^-ir r ^ C{u)ud-i{du). 

In particular iox C, = 1, W — the unit ball and d — 2 and o? = 3 we get the values 

32 

Vb2(Vo1i, Aj^o) = TT^ and Vsii^oh, Kso) = y^r^ 

respectively. 

Proof of Proposition 1 First, [9, Cor. 3.7] provides a formula for the second moment 
of the volume of the typical Poisson cell of a stationary Poisson hyperplane tessellation in 
W^. In terms of the zonoid 11 it reads 



:Vol^(TypicalCell(PHT(A))) = ^d^^^ 
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where we have used formula [40, (4.63)] . Moreover, the mean volume of TypicalCell(PHT(A)) 
is given by 

EVol,(IVpicalCell(PHT(A))) = ^^J-^ 

according to [40, Thm. 10.3.3 and (10.4)]. Using now Eq. (4.61) ibidem, key property (a) 
of STIT tessellations and replacing d by d — 1 in the last to formulas we obtain (72). □ 

5.2 Gaussian Limits in the Planar Case 

In the planar case the asymptotic behaviour of E^(F(1,H/r)) turns out to be Gaussian. 
We write 

t{s,R) := exp([logi?- log log i?](s - 1)) = R'-\logRy-' (73) 
and we define the total length process 

>Cf := -^f^,{Y{T{s, R), Wn)), s e [0, 1]. (74) 
rtVlog R 

The main result of this subsection is 

Theorem 6 The total length process {Cf'^)se[o,i] converges in law, as R ^ oo, on the 
space V[0,1] of cadlag functions on [0, 1] endowed with the usual Skorokhod topology, to 
(VVV(0,A)W,).e [0,1] where, again, W(.) stands for the standard Wiener process. 

In particular, for = Voli we have 

R^/^og R 

for the stationary and isotropic STIT tessellation 1^(1) in the plane, where =^ means 
convergence in law. 

Proof of Theorem 6 Note first that 

t{0,R) = t{1,R) = 1, j^T{s,R) ^T{s,R)[logR-loglogR]. (75) 

Thus, defining the auxiliary process 

^ = Ms \ =[f:^{Y{r{s, R), Wn)) - E^(r(T(0, R), Wn))] 

R^/\ogR-\og\ogR ^^^-^ 

and using (43) with Wr := RW and under variable substitution s := t{u,R) and t := s 
with LHS variables corresponding to the notation of (43) and RHS to that used here, we 
see that, by (75), 

{M,)U, and (m^ - f ^^A^.iY{r{u, R), Wn))du) (77) 
V Jo ^ / se[o,i] 
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are $5T-(s,ij)-martingales. In particular, see [16, Thm. 4.2], the predictable quadratic varia- 
tion process (M)^ is given by 

{M)s = J^^ ^^^A^.{Y{T{u,R),WR))du, s e [0,1]. (78) 

Repeating the argument leading to (68) we see that 

lim ^^^'^^ A^2(Y(t(s, R), Wr)) ^ Vw(4>, A) in probability, uniformly in s e [0, 1]. 

Note that the uniformity in s comes, as in the case of (68), from the relation (67) implying 
that, in distribution, all instances of the LHS for different values of s are just scaling 
instances of the same object ^"^^4^2(^(1, W^^)) for R — R/t{s,R) and thus, in terms of 
the considered convergence in probability to a deterministic limit, we are just dealing with 
a single asymptotic statement. Consequently, by (79) and in full analogy to (70), 

lim (M)^ = f Vw{(pA)du = sVw{(t>A) in probability (80) 
Jo 

Thus, we are again in a position to apply [48, Thm 2.1] yielding the functional convergence 

in law, as i? — )■ cxD, in I^[0, 1] of {Ms)s^[o,i] to {\/Vw{4>, ^)y^s)s&[o,i]- Indeed, condition 
(ii.6) there is just (80), condition (ii.4) is trivial in view of (78), whereas the condition 
(ii.5) is verified by noting that, with probability one, J{M] 1) = ^^^^^^ 0{R(\\a,m{W)) = 
0(l/\/logi?) so that in particular limfi^oo^>J^{M; 1) = as required. Denoting now by 
C^'^ the correction term (i?VIogTR)-^E0(F(r(O, i?), !¥«)) such that 



y log-R 

noting that logi? — log log ~ logi? and that, by the scaling property (c) of STIT tessel- 
lations and by (55), 

Var(C^'«) = 0{[R-'{\ogR)-'][Ry{\ogRf][{\ogRy{\og\ogR)]) = 0(loglogit:/logi?) 

(81) 

we see that the processes Mg and are asymptotically equivalent in ©[0, 1] as — >■ oo. 
This completes the proof of Theorem 6. □ 

Remcirk 6 In the context of proof of Theorem 6 it should he remarked that the 'negligible 

correction term' C^'^ has its variance of order O (log log i?/ logi?) and thus indeed tending 
to 0, hut extremely slowly. Consequently, although the Gaussian CLT holds for Cq'^ , 
it is quite natural to expect that the convergence rates are extremely slow, conjecturedly 
logarithmic. This is due to the fact that dimension 2 is the largest dimension (critical 
dimension) where the Gaussian limits are still present. In dimensions 3 and higher there 
is no Gaussian GLT and the 'correction term' analogous to C^'^ will turn out order- 
determining rather than negligible, as shown in the next subsection. 
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5.3 Non-Gaussian Limits for d > 2 



We claim that the argument of Subsection 5.2 above cannot be repeated for d > 2. Intu- 
itively, this is due to the fact that for d > 2 the variance order of T,^{Y{1, Wr)) is 0{R'^^'^~^^) 
and so even the very first faces born in the cell division process already do bring a non- 
negligible contribution to the overall variance. Thus, we cannot split the whole STIT 
construction into the warm-up phase {t G [0,R~^logR] for d = 2) with negligible variance 
contribution and the proper phase unfolding already in a typical STIT environment. In 
fact, we claim CLT does not hold for STIT length functionals in dimension greater than 
2! To see it, observe first that, by the scahng property (c) of STIT tessellations. 



R-^''-'%iY{l, Wr)) = W)). (82) 

Further, recall that by (43) the process R k-)- Ti^(Y{R,W)) is a square-integrable martin- 
gale with absolutely continuous predictable quadratic variation process given in (44) and, 
moreover, by (47) we have 

Ef:l{Y{R,W)) ^ r e{n{H))[ [ 1 " exp(-i?A(N])) ^^^^^^^^^ 
J[W] JwnH JwnH ^m^y\) 

which is bounded uniformly in R. Consequently, by the martingale convergence theorem, 
there exists a centered square-integrable random variable 5(VF) such that 

E{W) = lim E^(y(i?, W)) (83) 

a.s. and in and, moreover, 

VaTE{W)= f e{n{H)) f f —}—-dxdyK{dH). (84) 
J[W] Jw^^H Jw^^H ^^\[^y\) 

Using now (82) we readily conclude that 

i?-('^-i)E^(r(l, W)) =^ E{W) (85) 

as i? — )■ cxD. 

We claim that the variable S(H^) is not normal. Even though we are able to show this fact 
for all W and translation invariant A by establishing non-Gaussian tail decay, for simplicity 
we only give a proof for an easily tractable particular case, postponing the study of more 
involved properties of the random field S(VF), 1^ C R*^ for a future paper. Namely, take 
W = [0, 1]'' and 

^ ^= E / ^re,+4dr (86) 
i=l 

where e^, i — 1, . . . ,d are vectors of the standard orthonormal basis for R*^ and (^^e +e+ 
is the unit mass concentrated on the hyperplane orthogonal to Cj in distance r from the 
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origin. Consider the event £n, N > 0, that only hyperplanes orthogonal to ei have been 
born during the time [0, 1] of the MNW-construction and their number exceeds N. Observe 
that, in view of the form (86) of A, P(£jv) — exp(— c?) Yl'k=N+i h t^^s 

\og¥(SN)^-e(N\ogN), (87) 

where by ©(■) we mean something bounded both from below and above by multiplicities 
of the argument, i.e. 0(-) = O(-) fl ^{■)- Further, given a fixed collection of all hyper- 
planes Hi, . . . ,Hk, k > N, born at times between and 1, on the event we see that 
the conditional law of S(iy) coincides with that oi k — d plus sum of independent copies 
^1, . . . , $,k+i of S(W^i)) • • • ) S(M4'+i) respectively, where Wj, j = 1, . . . ,k + l, are paral- 
lelepipeds into which W is partitioned by i^i, . . . , H^.. Note that the extra k above is the 
sum of {d — l)-volumes of Hi (IW whereas —d = —'K'E^{Y{1, W)) is the centering term. 
Since Var[^i -|- . . . -|- Ck+i] = Var S(lV"j) which is bounded above by Var E{W) in view 

of (84), by Chebyshev's inequahty we get P(Ci + . . . > -2^VarS(l^)) > 3/4. Thus, 

in view of (87) 

F{E{W) >N)>If (V2V^^+.) = exp(-e(AriogAr)). (88) 

Since Gaussian variables exhibit tail decay of the order exp(— 0(A^^)), the random variable 
S(Vr) cannot be normal, which completes our argument. 

6 Comparison with other Tessellation Models 

This section is devoted to a comparison of the stationary and isotropic STIT tessellations 
Y{t) with stationary and isotropic Poisson hyperplane tessellations and Poisson-Voronoi 
tessellation with the same surface intensity t > stemming from a stationary Poisson pro- 
cess (we restrict ourselves to the isotropic case, since Poisson-Voronoi tessellations based on 
stationary Poisson processes are automatically isotropic). This is of particular interest for 
the statistical analysis of random tessellations, for example when constructing asymptotic 
confidence intervals for mean value estimators and for goodness-of-fit tests. A first example 
was considered in in [34] where the linear contact distribution was used as a criterion. The 
mean values obtained in [32], [33] and [46] can be used for such a test and now also some 
second order quantities. More details on goodness-of-fit tests for random tessellations can 
for example be found in [39] together with a couple of examples. 

Denote by PHT(t, W) the stationary and isotropic Poisson hyperplane tessellation inside 
a bounded convex window W C M.'^. For d = 2 we will write PLT{t,W) instead of 
PHT(t, ly), since hyperplanes are just lines in this case. It was shown in [12] that the 
asymptotic variance of the total /c-volume of the /c-faces of PHT(t, W) equals 

Yai{Vo\k{PRT{t^V)))^t^^''->'^-\d-ky(f]'4 f^V^' JiW) 
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with 

J{W)^ f Voll_,{WnH)Aiso{dH). (89) 
Jh 

In the case W = Bj^ this integral can be evaluated explicitly and we have 

J{B^) - (2d-l)! ^^^^ ■ 

For Poisson-Voronoi tessellations PVT(7, W) (here restricted to W) constructed from a 
stationary Poisson point processes in with intensity 7 we have from [13] 



Var(Volik(PVT(7, W))) - T^'^^7'^ ^olaiW). 

lyt 

Brakke [4] has obtained 



Unfortunately, no analj^ic expression for r^f^ is currently known. Only for the planar case 



rf ^ ?s 1.0445685 

as the result of the evaluation of a rather involved multiple integral. That is the reason, 
why we will restrict in this context to the planar case d — 2. Here we have 

J{Bl) = ^i?^ 

and 

1 a. 

Var(Voli(PLT(t,S|))) ~ —tR^. 

o 

For the planar Poisson-Voronoi tessellation PVT(7) it is well known that its edge length 
intensity La is related to 7 by = 2^/7. Taking now La = t a& the parameter, we obtain 

Var(Voli(PVT(iV4,S^)) ~ 7^Tf^i?^ 

which is in particular independent of the edge length intensity t. The last fact can in- 
tuitively explained as follows: per unit area we have, in mean, ^ cells with mean total 
perimeter ci • t. Thus, a single cell has mean perimeter c^t and Variance c^t^. Since the 
cells are 'almost independent' of each other (sec [40, Chap. 10.5]), these variances add up 
and cancels out (here Ci, C2, C3 stand for some real constants). In fact a lot more is known 
about the typical geometry of Poisson-Voronoi tessellations in the plane, including further 
explicit information, see e.g. Calka [6, 7]. In Fig. 3 the graphs of the total variance of the 
edge length of Y{1, S|), PLT(1, B\) and PVT(l/4, B\) as a function of R are compared 
and the behaviour discussed above is visualized. As already observed in earlier papers, 
STIT tessellations turn out to interpolate between Poisson line and Poisson-Voronoi tes- 
sellations in some sense. The list where this is true can now be extended by the variance 
of the total edge length (or more generally the total surface area) . 
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PLT PVT STIT 




Figure 3: Comparison of the asymptotic variance of the total edge length for a STIT 
tessellation Y{1) (thick hne) with that of a Poisson hne tessellation (dashed line) and a 
Poisson-Voronoi tessellation (thin line) with the same surface intensity in a growing ball 



We can also compare the reduced second moment function and the pair-correlation func- 
tion of Y{t) given by (61) and (62) with the corresponding functions of a stationary and 
isotropic Poisson hyperplane tessellation PHT(i) with the same surface intensity t > 0. 
Wc will denote these two functions by K^^'^^^\r) and g^^^^^\r). From Slivnyak's theorem 
for Poisson processes one immediately infers that tK^^^^^\r) is given by tK^^'^^^\r) — 
Yo\d-i{B^-^) + tYo[d{B^) and thus 

From the definition of the pair-correlation function we obtain now 

^ ^ dKar'i-^ dr dK^tr ' 

Especially for the planar case d — 2, i.e. for the Poisson line tessellation PLT(t) we have 
the formulas 

An ^'{r) — T:r-\ and On ^'{r) — l^ . 

We have to skip the comparison of the K-function and the pair-correlation function be- 
tween the STIT and Poisson-Voronoi tessellations, since nothing seems to be known about 
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Figure 4: Fig. 3: Comparison of the i^'-function (left) and the pair-correlation function 
(right) for a STIT tessellation Y{1) and a Poisson Une tessellation PLT(l) with the same 
parameter 



the pair-correlation function of the random length measure of a planar Poisson- Voronoi 
tessellation. However, the second order analysis of the point process of nodes of a Voronoi 
tessellations is rather involved and we refer the reader to Heinrich and Muche [14] and 
references therein for details, including numerics. In Fig. 4 we can see that the K-function 
and the pair-correlation function of a planar STIT tessellation and a Poisson line tes- 
sellation are very close together and it seems - having statistical applications in mind - 
that comparing variances could be more fruitful and could yield much better information 
about the underlying tessellation model than a comparison of the reduced second moment 
function or the pair correlation function (the picture is essentially the same for any space 
dimension) . 

In the last 15 years, central limit theorems for random tessellations have been considered 
by many authors, mostly for the case of stationary and isotropic Poisson hyperplane or 
stationary and isotropic Poisson- Voronoi tessellations, see Heinrich [12] and the references 
therein. Elegant proofs of central limit theorems for Poisson hyperplane tessellations rely 
on the fact, that many functional may be expressed as U-statistics with a special kernel 
function. Apphcation of Hoeffings's decomposition and Hoeffding's central limit theorem 
for U-statistics leads now to hmit theorems for geometric functionals of Poisson hyperplane 
tessellations. The classical Cramer- Wold device also allows multivariate extensions. Note 
that in all cases the limit distribution is a normal distribution. Note that the non- additive 
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ovoid functional J{W) from (89) may also be written as 



where Id is the d-th chord power integral of W (here one uses again the affine Blaschke- 
Petkantschin formula to show the equality) . It means that the asymptotic variance of the 
total surface area of a Poisson hyperplane tessellation depends on W by Id{W) and also 
on the surface intensity t. 

The central limit theorems for functionals of the Poisson- Voronoi tessellations, as provided 
by many authors [1, 12, 13, 14, 35, 36, 37, 38], are mainly based on their strong mixing 
properties and the fact that they exhibit exponential decay of dependencies, whereas Pois- 
son hyperplane tessellations are strongly mixing, but long-range dependent, which means 
that we observe slow decay of the correlations between distant parts of the tessellations. 
Unfortunately - as already explained above - the results are much less explicit in this case, 
since most often analytic expressions for variances are not known. But again, the limit 
distribution were shown to be normal distributions in all cases. The asymptotic variances 
are shown to depend on the shape of the sequence of growing observation windows Wr only 
through their volume Yold{W) which may itself also be (artificially) expressed as -^Ii{W), 
where Ii{W) is the first-order chord power integral of W. They also depend (beside the 
planar case) on the parameter t. 

In the following table the facts about asymptotic variances of the total surface area of a 
PHT, a PVT and a STIT tessellation in the sequence of growing convex observation win- 
dows Wr C are summarized for the isotropic cases (c^ denotes some space dimension 
dependent constant, but not necessarily always the same, and Ik{W) is the k-th chord 
power integral of W). 



Model 


PVT 


STIT 


PHT 


Var(Volrf_i(-,W^fl)) 


Cdt'^-'^h{W)R'^ 


C2h{W)RHogR {d = 2) 


CdtId{W)R'd-^ 



There are also some central limit theorems for stationary iterated tessellations, which were 
motivated by potential applications in telecommunication. However, these theorems do not 
lead to corresponding results for STIT tessellations. As a nice survey for the central limit 
theory of Poisson hyperplane, Poisson- Voronoi and iterated tessellations with full proofs 
and several applications we recommend the small book [39]. Another recommendable 
source is the recent survey [8] . 

References 

[1] Baryshnikov, Y.; Yukich, J.E.: Gaussian limits for random measures in geometric 
probability, Ann. Appl. Probab. 15, 213-253 (2005). 



48 



[2] BiLLlNGSLEY, P.: Convergence of probability measures, Wiley (1968). 

[3] BOGACHEV, V.I.: Measure Theory, Volume 1, Springer (2007). 

[4] Brakke, K.A.: Statistics of random plane Voronoi tessellations, unpublished manuscript 
(1987). 

[5] Breiman, L.: Probability, SIAM (1992). 

[6] Calka, p. Precise formulae for the distributions of the principal geometric characteristics 
of the typical cells of a two-dimensional Poisson-Voronoi tessellation and a Poisson line 
process. Adv. Appl. Probab. 35, 551-562 (2003). 

[7] Calka, P. An explicit expression for the distribution of the number of sides of the typical 
Poisson-Voronoi cell. Adv. Appl. Probab., 35, 863-870 (2003) 

[8] Calka, P. Tessellations, In New Perspectives in Stochastic Geometry, edited by W.S. 
Kendall and I. Molchanov, Oxford University Press (2009). 

[9] Favis, W.; Weiss, V.: Mean values of weighted cells of stationary Poisson hyperplane 
tessellations in W^, Math. Nachr. 193, 37-48 (1998). 

[10] Georgii, H.-O.: Gibbs measures and phase transitions, de Gruyter (1988). 

[11] Harris, J.W.; Stocker, H.: Handbook of Mathematics and Computational Science, 
Springer, New York (1998). 

[12] Heinrich, L.: Central limit theorems for motion-invariant Poisson hyperplanes in expand- 
ing convex windows. Preprint (2008). 

[13] Heinrich, L.; Schmidt, H.; Schmidt, V.: Limit theorems for functional on the facets 
of stationary random tessellations, Bernoulli 13, 863-891 (2007). 

[14] Heinrich, L.; Muche, L.: Second-order properties of the point process of nodes in a 
stationary Voronoi tessellation. Math. Nachr., 281, 350-375 (2008) 

[15] Itoh, H.: An analytical expression of the intercept length distribution of cubic particles. 
Metallography 3, 407-417 (1970). 

[16] Jacod, J., Shiryaev, A.N.: Limit theorems for stochastic processes, Grundlehren der 
Mathematischen Wissenschaften 288, 2nd ed.. Springer (2003). 

[17] Kato, T.: Perturbation Theory for Linear Operators, Springer (1966). 

[18] KiPNis, C, Landim, C: Scaling Limits of Interacting Particle Systems, Springer (1999) 

[19] Karatzas, I., Shreve, S.E.: Brownian motion and stochastic calculus. Graduate Texts 
in Mathematics, Second Ed., Springer (1998). 

[20] Liggett, T.M.: Interacting particle systems, Springer (1985). 



49 



[21] Lachieze-Rey, R.: Strong mixing property for STIT tessellations, arXiv:0905.1145v3 
[math.PRj (2009). 

[22] Matheron, G.: Random Sets and Integral Geometry, Wiley (1975). 

[23] Mattila, p.: Geometry of sets and measures in Euclidean spaces: fractals and rectifiability, 
Cambridge University Press (1995). 

[24] Mecke, J.: Joint distribution of direction and length of the typical I -segment in a homo- 
geneous random planar tessellation stable under iteration. Journal of Contemporary Math- 
ematical Analysis 44, 45-53 (2009). 

[25] Mecke, J.; Nagel, W.; Weiss, V.: Length distributions of edges in planar stationary 
and isotropic STIT tessellations, Izvestija Akademii Nauk Armenii, Matematika, 42, 39-60 
(2007). 

[26] Mecke, J.; Nagel, W.; Weiss, V.: A global construction of homogeneous random planar 
tessellations that are stable under iteration, Stochastics 80, 51-67 (2008). 

[27] Mecke, J.; Nagel, W.; Weiss, V.: The iteration of random tessellations and a construc- 
tion of a homogeneous process of cell division. Adv. Appl. Prob. 40, 49-59 (2008). 

[28] Miles, R.E.: Poisson flats in Euclidean spaces. Part II: Homogeneous Poisson fiats and 
the complementary theorem. Adv. Appl. Probab. 3,1-43 (1971). 

[29] Miles, R.E.; Mackisack, M.S.: A large class of random tessellations with the classic 
Poisson polygon distribution. Forma 17, 1-17 (2002). 

[30] Nagel, W.; Weiss, V.: Limits of sequences of stationary planar tessellations, Adv. Appl. 
Prob 35, 123-138 (2003). 

[31] Nagel, W.; Weiss, V.: Crack STIT tessellations: characterization of stationary random 
tessellations stable with respect to iteration. Adv. Appl. Prob 37, 859-883 (2005). 

[32] Nagel, W.; Weiss, V.: STIT tessellations in the plane, Rendiconti del circulo matematico 
di Palermo, Serie II, Suppl. 77, 441-458 (2006). 

[33] Nagel, W.; Weiss, V.: Mean values for homogeneous STIT tessellations in 3D, Image 
Anal. Stereol. 27, 29-37 (2008). 

[34] Nagel, W.; Mecke, J.; Ohser, J.; Weiss, V.: A tessellation model for crack patterns 
on surfaces. Image Anal. Stereol. 27, 73-78 (2008). 

[35] Penrose, M.D.: Gaussian limits for random geometric measures. Electronic Journal of 
Probability, 12, 989-1035 (2007) 

[36] Penrose, M.D.: Laws of large numbers in stochastic geometry with statistical applications, 
Bernoulli, 13, 1124-1150 (2007) 



50 



[37] Penrose, M.D.; Yukich, J.E.: Central limit theorems for some graphs in computational 
geometry, Ann. Appl. Probab., 11, 1005-1041 (2001) 

[38] Penrose, M.D.; Yukich, J.E.: Weak laws of large numbers in geometric probability, Ann. 
Appl. Probab., 13, 277-303 (2003). 

[39] Schmidt, H.: Asymptotic Analysis of Stationary Random Tessellations, VDM Verlag Dr. 
Miiller, Saarbrucken (2008). 

[40] Schneider, R.; Weil, W.: Stochastic and Integral Geometry, Springer, Berlin (2008). 

[41] Schneider, R.; Weil, W.: Classical stochastic geometry. In New Perspectives in Stochastic 
Geometry, edited by W.S. Kendall and I. Molchanov, Oxford University Press (2009). 

[42] Stoyan, D.: Inequalities and bounds for variances of point processes and fibre processes. 
Math. Operationsf. Statist., Ser. Statistics 14, 409-419 (1983). 

[43] Stoyan, D.; Kendall, W.S.; Mecke, J.: Stochastic Geometry and its Applications, 
Second Edition, Wiley, Chichester (1995). 

[44] Stoyan, D.; Stoyan H.: Fractals, Random Shapes and Point Fields, Wiley, Chichester 
(1994). 

[45] Thale, C, Moments of the length of line segments in homogeneous planar STIT tessella- 
tions. Image Anal. Stereol. 28, 69-76 (2009). 

[46] Thale C; Weiss, V.: New mean values for spatial homogeneous random tessellations 
stable under iteration, in preparation. 

[47] Weiss, V.; Ohser, J.; Nagel, W.: Second moment measure and K -function for planar 
STIT tessellations submitted for publication. 

[48] Whitt, W.: Proofs of the maHingale FCLT, Probability Surveys, 4, 268-302 (2007). 



51 



